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Measurements of the penetration tensions in paper 
of a series of creosotes distilled from a single coal tar are 
reported (the penetration tension being the product of the 
surface tension and the cosine of the contact angle). It is 
found that these data conform to two equations: one, 
previously developed by others from theoretical considera- 
tions of the effects of variation in size and distribution of 
capillary pores upon the theoretical equation of flow, the 
other, an empirical relationship which is less difficult to 
apply in more practical problems where many data must 


be analyzed. The creosote-paper contact angles are shown 
to be zero or extremely small. Possible errors in the paper 
strip method are discussed and the importance of calibrat- 
ing each strip with a reference liquid is demonstrated. 
The failure of Washburn’s equation, which has received 
considerable attention in the field of fibrous penetration, 
is shown to be due to the neglect of the hydrostatic head 
term. The penetration tensions of these creosotes are 
unaltered by the presence of adsorbed water. 





INTRODUCTION 


HILE industry has long been fully cog- 

nizant of the commercial importance of 
solid-liquid relations in the various processes of 
manufacture which depend upon the movement 
of liquid penetrants in fibrous materials, the laws 
governing the dynamics of liquid flow in such 
materials are still in the early stages of develop- 
ment. This has not been due to a lack of research 
on the subject but rather to the inherent difficul- 
ties encountered in a study of this nature. As 
examples of such difficulties may be cited the 
inability to duplicate samples; variations in the 
size, shape and distribution of the capillary 
pores; moisture effects; possible swelling brought 
about by the adsorption of the liquid penetrant; 
visual difficulties and a lack of knowledge of the 
specific mutual relations at solid-liquid inter- 
faces. As a consequence, the equations of flow 
in fibrous materials have more generally been 
derived from a consideration of liquid flow in 
uniform capillaries. 
~ * This paper is one of a series entitled ‘‘ Chemical Studies 


of Wood Preservation” being published in this journal 
and in Industrial and Engineering Chemistry. 


The most direct method which has been em- 
ployed to study the penetration of liquids in 
nonuniform capillary systems is known as the 
Klemm test. It depends upon the fact that 
wetting liquids, due to capillarity, will rise in 
strips of a fibrous material when one end of the 
strip is brought in contact with a liquid. 

Impregnation of wood is the fundamental 
concern of this investigation but the difficulties 
involved in the preparation and study of thin 
uniform strips of wood suggested the substitu- 
tion of a more convenient fibrous material. The 
somewhat similar composition and structure that 
exists in paper and wood encouraged the belief 
that useful information could be obtained by a 
study of the penetration of creosote into paper 
strips. 

The eight creosotes of the Fulweiler series were 
selected as penetrants since the necessary phys- 
ical constants had already been determined.'! 
These creosotes were all distilled from a common 
tar but differed in volatility because of variations 
in the temperature of distillation. 


1 Frosch, Physics 6, 165 (1935). 
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HISTORICAL 


The dynamic movement of liquids in capillary 
systems such as soil and blotting paper was con- 
sidered by Bell and Cameron.? They showed 
that the penetration could be expressed by the 
empirical relation. 


h" =kt, (1) 


in which ¢ is the time required to penetrate the 
structure to a height 4 and 7 and k are constants. 
They found the value of m to be 2.3. 

Later Ostwald* deduced the equivalent equa- 
tion 


h=kt", (2) 


which Lucas‘ applied to the rise of various liquids 
in filter paper and found values of m ranging from 
0.424 to 0.500 corresponding to exponents of 
2.00 to 2.36 in Eq. (1). Simmonds arrived at a 
value of m equal to 0.425 for the penetration of 
water into blotting paper, cardboard and un- 
beaten sulphite pulp. 

More recently Washburn® considered the 
penetration of liquids in uniform capillaries. 
He showed that inertia effects are inappreciable 
after the penetration has progressed a measur- 
able distance and, for small capillaries, may be 
neglected. He also considered hydrostatic effects 
but suggested that these may be neglected in 
comparison with capillary forces when dealing 
with small capillaries. His equation then takes 
the form 


h?=rty cos 0/2n; (3) 


neglecting the viscosity of the displaced gas, 
where y is the surface tension of the liquid, 7 is 
the viscosity, 7 the radius of the capillary and @ 
is the angle of contact between the liquid and the 
solid. Whitehead and Greenfield’ applied this 
equation to the capillary rise of oils in cable 
insulation papers and found that a particular 
value of r for each paper satisfies the penetration 
data. 


? Bell and Cameron, J. Phys. Chem. 10, 658 (1906). 
* Ostwald, Kolloid, Zeits. 2, Suppl. No. 11 (1908). 

* Lucas, Kolloid Zeits. 23, 15 (1908). 

* Simmonds, Paper Trade J. 97, 40 (1933). 

* Washburn, Phys. Rev. 17, 273 (1921). 

* Whitehead and Greenfield, Physics 3, 324 (1932). 
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Peek and McLean have shown that while the 
effect of hydrostatic head may be small, it is not 
always insignificant even for the penetration of 
liquids into small capillaries. They considered 
not only the effects of hydrostatic head in their 
study but also the effects of variations in size and 
distribution of pores upon the theoretical equa- 
tion of flow for uniform capillaries. They de- 
rived an equation in the form 


dh Aycos@ Bdg 


’ 4) 
dt 4nh 8n ( 


where A and B are constants for each paper strip 
depending upon the pore size distribution func- 
tion. They have shown that B will be equal to 
A? only for the case of uniform capillary pores 
throughout the paper structure when A becomes 
analogous to r in Eq. (3). The application of Eq. 
(4) to experimental data is by graphic means. 


FIBROUS PENETRATION 


The kinetics of flow involved in the capillary 
rise of liquids in fibrous structures have already 
been considered by a great many investigators 
but will be briefly discussed here for the sake of 
completeness. 

When acapillary of uniform radius 7 is brought 
in contact with a liquid which will wet it, the 
liquid will rise in the tube because of capillary 
pressure. This motion will be opposed by the 
hydrostatic head of the liquid so that the total 
force per unit cross section producing motion 
will be 

2y cos 6 r—hdg (5) 


for a vertical capillary. According to Poiseuille’s 
law, the rate at which a liquid will penetrate a 
capillary, neglecting the viscosity of the dis- 
placed gas, is equal to 


Pr 8nh, (6) 


where P is the total force producing motion. 
Substitution of the capillary tension and hydro- 
static head forces for the value of P leads to a 
relationship very similar to Eq. (4). When the 
rate of flow is equal to zero, the equation has the 


* Peek and McLean, Ind. Eng. Chem., Anal. Ed. 6, 255 
(1934). 
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form of Jurin’s law which is used in a static 
method for the measurement of surface tension. 
The differential form of Eq. (3) results when the 
hydrostatic head of the liquid is neglected. 

The penetration of a liquid in a paper strip is 
obviously more complex than that in uniform 
capillaries. The paths of liquid flow have become 
a maze of branching and connecting capillaries 
varying greatly in size throughout their length. 
Peek and McLean® considered the effects of 
variations in size and distribution of capillary 
pores upon the rate of penetration of liquids in 
paper strips. They assumed that the same distri- 
bution function of pore sizes holds throughout 
the length of the paper strips. Such an assump- 
tion is indeed reasonable when it is considered 
that the rate of liquid penetration in a given 
paper strip is substantially the same regardless 
of which end of the strip is submerged beneath 
the liquid. This obviously could not be the case 
unless the various pore sizes were distributed 
with substantially the same frequency through- 
out the paper. A constant distribution of pore 
sizes throughout the length of a paper strip 
requires a constant pore volume in any cross 
section of the paper. This seems to be indicated 
by the uniformity in density of small cross- 
sectional elements of a given paper strip. 

It was in accordance with the above reasoning 
that Peek and McLean’ deduced Eq. (4). They 
found that their equation gave a very good fit 
to their penetration data. They also found that 
A and B did not vary except for experimental 
error among all of the paper strips which they 
studied. It will be shown subsequently that A 
and B vary considerably among the paper strips 
considered herein but are essentially constant 
for a given paper strip when different liquids are 
employed. Since A cos 6 and B can be determined 
from a single penetration study, the problem 
resolves itself into the determination of A 
for each paper strip either by an independent 
method or by the penetration of a reference 
liquid. 

The very nature of the structure in a paper 
strip would necessitate the determination of A 
by penetration studies. Efflux measurements 
utilizing either gases or liquids are impossible as 
the capillaries in the direction to be measured 
are parallel to the length of the paper. The 


capillary penetration of a reference liquid can 
be used to determine A if the contact angle be- 
tween the paper and liquid is known. 

Other investigators have shown that in general 
the contact angle depends upon the liquid used 
and the surface condition of the solid. The direct 
measurement of paper-liquid contact angles in 
paper capillaries obviously cannot be made. 
However, it should be possible to find a liquid 
whose rate of penetration would not vary with a 
change in the surface condition of the paper 
capillaries thus providing evidence of a zero 
contact angle in all instances. The conditioning 
of the surface must be accomplished by a method 
which will not alter the pore size by chemical or 
mechanical means and thus invalidate the re- 
sults. McLean and Kohman® have shown that 
oxidation of paper does not affect the rate of 
penetration of benzene whereas it decreases the 
rate of penetration of water. This demonstrates 
that the paper-water contact angle has probably 
been changed; while that of benzene has re- 
mained unaltered, for oxidation would hardly be 
expected to change the structure of the paper, at 
least in such a way as to compensate exactly for 
a change in contact angle. Unpublished data 
have also shown the rate of penetration of 
benzene to be unaltered with other types of 
surface conditioning. A constant angle of contact 
with a changing surface condition would prob- 
ably be equal to zero, but even if this were not 
the case, the assumption of a zero contact angle 
would still give results of practical importance. 
The constant error introduced in the A values 
for each paper strip by assuming a zero contact 
angle for the reference liquid when it is other than 
zero would alter the absolute values of the 
penetration constants of other liquids but not 
their relative values. 

In view of these considerations, when a paper 
strip is calibrated with benzene and a second 
liquid is measured over the calibrated portion 
of the strip, the penetration constants relative to 
benzene can be determined by Eq. (4). It is 
obvious that the above equation cannot be ex- 
pected to apply if the effective paths of flow are 
different for benzene and the second liquid. 
This condition may be brought about by swelling 


® McLean and Kohman, Elect. Eng. 53, 255 (1934). 
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of the fibrous structure or by mechanical dis- 
persion of the loosely bound paper fibers. There 
seems to be no available information in the 
literature concerning the change in pore size of 
paper during impregnation with various liquids 
although Roth!® has shown that the cell walls of 
pine wood increase in thickness with adsorption 
of water, i.e., the pore size diminishes. This has 
more recently been demonstrated by Stamm" 
who measured the flow of air at various relative 
humidities through wood sections. Since wood is a 
continuous structure and is therefore not sub- 
ject to mechanical dispersion by liquid pene- 
trants, each value of r in a general formulation 
must be replaced by (r—or) where o is the frac- 
tional decrease in the radius of each capillary 
brought about by swelling and would be equiva- 
lent to the internal swelling coefficient. It can be 
safely assumed that benzene does not swell pine 
appreciably, for Hasselblatt’® found that benzene 
swells birchwood less than one percent externally, 
a figure which is believed to be approximately 
the same for internal swelling. 

Paper, unlike wood, is a discontinuous struc- 
ture of innumerable cellulosic fibers distributed 
in a more or less random manner and being rather 
loosely bound together. If these cellulosic fibers 
are assumed to be round incompressible rods, 
they might be arranged either in a closely packed 
incompressible condition or in a loosely packed 
compressible structure. The former would ob- 
viously result in increased pore size with swelling 
of the rods whereas the latter would more likely 
result in decreased pore size. It is believed that 
the loosely packed structure probably occurs in 
paper since pressure treatments followed by 
calendering are known to reduce the porosity of 
paper. It would appear from this that the effect 
of swelling liquids on paper would be similar to 
that on wood; but this is not necessarily the case, 
for mechanical dispersion of such a loosely 
bound structure may and probably does obtain 
during impregnation with liquids. The effects of 
mechanical dispersion may be larger than those 
of swelling in low density papers but may be 
less in more dense papers. Whether the pore size 


© Roth, Bot. Ztg. 52, 217 (1894). 

"' Stamm, Physics 6, 334 (1935). 
' 8 a Zeits. f. anorg. allgem. Chemie 154, 375 
1926). 
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is increased or decreased depends therefore upon 
the type of paper and liquid used. It is also 
evident that with certain types of papers both 
effects might be alternately predominant which 
might cause fluctuations in the values of the 
penetration constants when check runs are 
attempted. Such fluctuations have been noted 
when dealing with water and other liquids which 
are known to swell cellulose but they are not 
reported here since the scope of the present 
investigation does not deal with the changes in 
pore size brought about by the penetration of 
liquids in paper strips. The possible effects of 
these have been considered merely to point out 
one of the inherent sources of error in the paper 
strip method. 


EXPERIMENTAL 


The creosotes were held in an open dish in the 
original determinations by the paper strip 
method but it was found that, due to evaporation 
of volatile constituents, the viscosities had 
changed by as much as three percent after a 
single determination of 40°C. These results were 
disregarded and an apparatus designed which 
would remain closed throughout the course of a 
determination. This apparatus is shown in Fig. 1 
with a paper strip in position for a penetration 
study. The paper is fastened to hooks on a 
glass rod by means of paper clips soldered to wire 
eyelets. A spring is added at the upper end to 
facilitate mounting and to keep the strips under 
a slight tension. The glass rod is suspended by a 
thread from a glass windlass by means of which 
it can be raised or lowered. Guides arrest rota- 
tion of the mounting so that a lamp placed 
behind the apparatus will furnish a clear view of 
the liquid meniscus at all times. The system can 
be evacuated by applying suction to the venting 
stopcock. Constant moisture conditions are 
maintained within the system by placing hu- 
midity controls such as dehydrite or saturated 
sodium chloride solution in the side bulb. By 
means of this arrangement the paper strip can be 
brought in contact with the liquid to be meas- 
ured and the measurements made _ without 
disturbing the established equilibrium. The 
apparatus thus assembled is held vertically in 
an oil bath controlled at the desired temperature 














' \v —_ | eel ‘vy o6VvSs _- — ball Dall a = 


—-_— We fF 








REEL 








MERCURY 
" SEAL 
‘y 
DRYING AGENT it 
OR ia 
SALT SOLUTION Hl 
m 
see 
gs! 
| 
eos att} PAPER 
' STRIP 











LIQUID 








GUIDES 


Fic. 1. Apparatus for studying penetration of liquids into 
paper strips. 


with an accuracy of +0.05°C by means of a 
mercury column thermostatic control. 

In studying the rate of penetration into paper, 
strips approximately one cm wide are cut from a 
section of Whatman No. 4 filter paper. Each 
strip is lightly marked at half-cm intervals and 
successive markings are numbered. The strips 
are extracted with benzene to remove any 
extractives which might affect later results. 
Benzene is removed by drying at 60°C in vacuum 
and the samples are then conditioned in evacu- 
ated desiccators for at least one month. The 
humidities are controlled with dehydrite or 
saturated sodium chloride solution. 

A sufficient quantity of benzene is placed in 
the apparatus with the side bulb containing one 
of the humidity controls. A paper strip is 
quickly transferred from the desiccator to the 
apparatus and slight suction is applied to prevent 
disruption of the mercury seal by the pressure 
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developed with rise in temperature. The system 
is allowed to remain at constant temperature for 
at least 24 hours before the measurements are 
made. The paper is brought into contact with the 
benzene and the times of rise are taken at suc- 
cessive half-cm intervals. The paper is removed, 
dried at 60°C in vacuum and reconditioned in 
the original desiccator for another month. More 
than one determination with benzene was made 
on some of the papers to determine the repro- 
ducibility of the data. 

The strips thus calibrated are used in the above 
manner to determine the times of rise of the 
Fulweiler creosotes at various temperatures. 
The data for benzene are all at 40°C but it is be- 
lieved that the change in pore size due to ex- 
pansion with increasing temperature is negligible. 
Whitehead and Greenfield? showed that pore 
size increased markedly with temperature in- 
crease but they do not mention any humidity 
control. Their relative humidity would be ap- 
preciably less at the higher temperatures and 


TABLE I. Some physical properties of the Fulweiler creosotes. 

















Sur- TIME oF RIsE** 
REsI- FACE 
DUE* TEN- (Seconds) 
SamM- Pa- ABOVE VIs- SION 
PLE PER 355°C Temp. cosity (dynes/ BEN- CREO- 
No. No. (Percent) (°C) (poises) cm) ZENE  SOTE 
F-1 1 43.84 40 0.1064 38.35 174 2673 
2 60 0.0544 36.41 149(a) 1133 
3 80 0.0320 34.54 210 1036 
q 100 0.0218 33.66 194 615 
F-2 5 40.25 40 0.0973 37.70 210 3008 
F-3 6 31.11 40 0.0675 37.59 210 1917 
7 60 0.0377 35.63 193 1140 
8 80 0.0243 33.85 168(a) 670 
9 100 0.0170 32.24 190 526 
F-4 10 28.52 40 0.0667 37.46 194 1897 
F-5 11 20.85 40 0.0525 37.41 156 1184 
12 60 0.0310 35.57 196(a) 919 
13 80 0.0211 33.90 197 617 
14 100 0.0148 31.99 229(b) 560(b) 


F-6 15 15.49 40 0.0445 37.49 194 1318 
F-7 16 12.21 40 0.0396 37.52 191 1023 


17 60 0.0247 35.41 201 780 

18 80 0.0170 33.64 190 544 

19 100 0.0125 31.39 199 459 
F-8 20 8.49 40 0.0344 37.01 190 968 

21 40 0.0344 37.01 199 953 

22 40 0.0344 37.01 184(c)  934(c) 








* Standard distillation of creosote, Standard 1le, American Wood- 
Preservers’ Association Manual of Recommended Practice. 

** Time of rise between three and ten centimeters. 

(a) Average of two or more benzene calibrations on the same paper 
strip. 

(b) Average of two determinations on different paper strips. 

(c) These determinations were made in equilibrium with saturated 
sodium chloride solution whereas the remaining ones were made in 
equilibrium with dehydrite. 
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consequently the amount of water adsorbed by 
the paper at the higher temperatures would be 
substantially reduced. The decrease in the 
amount of adsorbed water at the higher tempera- 
tures and consequent decrease in swelling might 
account for the large increase in pore size with 
increasing temperature which they found. 

The added precaution is also taken, in measur- 
ing the rates of rise of the creosotes, to make the 
point of initial contact for the creosote coincide 
with that of the benzene calibration. This limits 
the penetration data for the creosotes to cali- 
brated areas of the paper strips. 

Paper strips were discarded after a single 
determination with any of the creosotes since 
long extraction of the papers with benzene still 
leaves the strips slightly discolored. This is 
probably due to large molecules insoluble in 
benzene and free carbon which are present in 
coal tar creosotes. 

The penetration data thus determined are far 
too numerous for inclusion herein. The times of 
rise of benzene and the creosotes from 3 to 10 cm, 
as given in Table I, indicate the general trend of 
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the data. The values of 3 and 10 cm are chosen 
as it is felt that initial inertia and final equilib. 
rium effects are unimportant between these 
limits. 

The surface tensions, viscosities and residues 
not distilling up to 355°C previously determined! 
for the Fulweiler creosotes are also included in 
Table I for reference purposes. 

Repeated calibration of the paper strips with 
benzene gave successively lower values for the 
time of rise to the same height although this 
decrease in time was lessened by each successive 
determination. The explanation for this in- 
creased rate of penetration may lie in a general 
loosening of the cellulosic fibers as was previously 
discussed although it seems strange that the long 
extraction with benzene originally made on the 
papers did not eliminate this effect. The increased 
rates of penetration for benzene were not ex- 
ceedingly large and only the average calibration 
values were included in Table I. The magnitude 
of these variations is clearly indicated in the 
slope values for benzene given in Table IT. I 

It had been hoped that the effective radii of 


TABLE II. Penetration constants according to Eq. (4). 

















Temp. Run SLOPE INTERCEPT Ay cos 8 BX 108 Temp. RuN_ SLOPE INTERCEPT Ay cos 0 Bx108 
Liquip (°C) No. (cm*/sec.) (cm/sec.) (dynes) (cm?) LiguIp (°C) No. (cm?/sec.) (cm/sec.) (dynes) (em?) 
Benzene 40 1 0.399 —0.0155 0.00785 72.56 Benzene 40 1 0.254 —0.0028 0.00500 13.11 
F-1 40 2 0.024 —0.0010 0.01021 79.78 F-5 80 2 0.079 — 0.0009 0.00667 14.98 
Benzene 40 1 0.352 —0.0095 0.00693 44.47 Benzene 40 1 0.279 —0.0145 0.00549 67.87 
Benzene 40 2 0.380 —0.0083 0.00748 38.85 F-5 100 2 0.101 —0.0028 0.00598 33.19 
Benzene 40 3 0.386 —0.0083 0.00760 38.85 
F-1 60 4 0.046 —0.0010 0.01001 41.25 Benzene 40 1 0.274 —0.0045 0.00539 21.06 
Benzene 40 2 0.285 —0.0048 0.00561 22.47 
Benzene 40 1 0.337 —0.0161 0.00663 75.36 Benzene 40 3 0.290 —0.0045 0.00571 21.06 
F-1 80 2 0.066 —0.0029 0.00845 71.38 F-5 100 4 0.103 —0.0029 0.00610 34.37 
Benzene 40 1 0.269 —0.0045 0.00529 21.06 Benzene 40 1 0.305 —0.0099 0.00600 46.34 
F-1 100 2 0.080 —0.0012 0.00698 20.42 F-6 40 2 0.047 —0.0017 0.00837 58.53 
Benzene 40 1 0.270 —0.0045 0.00531 21.06 Benzene 40 1 0.276 —0.0056 0.00543 26.21 
F-2 40 2 0.018 —0.0004 0.00701 29.29 F-7 40 2 0.048 —0.0007 0.00760 21.53 
Benzene 40 1 0.237 —0.0033 0.00466 15.45 Benzene 40 1 0.246 —0.0033 0.00484 15.45 
F-3 40 2 0.025 —0.0003 0.00675 15.43 F-7 60 2 0.066 —0.0010 0.00652 19.45 
Benzene 40 1 0.264 —0.0040 0.00520 18.72 Benzene 40 1 0.286 —0.0064 0.00563 29.96 
F-3 60 2 0.046 —0.0008 0.00694 23.29 F-7 80 2 0.101 —0.0021 0.00687 28.56 
Benzene 40 1 0.296 —0.0050 0.00583 23.41 Benzene 40 1 0.281 —0.0059 0.00553 27.62 
Benzene 40 2 0.316 —0.0055 0.00622 25.75 F-7 100 2 0.117 —0.0023 0.00585 23.35 
F-3 80 3 0.075 —0.0011 0.00729 20.92 
Benzene 40 1 0.243 —0.0024 0.00478 11.23 
Benzene 40 1 0.275 —0.0051 0.00541 23.87 F-8 40 2 0.049 —0.0003 0.00674 8.09 
F-3 100 2 0.097 —0.0015 0.00660 20.25 
Benzene 40 1 0.256 —0.0024 0.00504 11.23 
Benzene 40 1 0.293 —0.0079 0.00577 36.98 F-8 40 2 0.051 —0.0005 0.00702 13.48 
F-4 40 2 0.030 —0.0007 0.00800 35.64 
Benzene 40 1(a) 0.286 —0.0055 0.00563 25.75 
Benzene 40 1 0.329 —0.0061 0.00647 28.55 F-8 40 2(a) 0.054 —0.0008 0.00743 21.56 
F-5 40 2 0.042 —0.0007 0.00882 28.18 
Benzene 40 1 0.254 —0.0054 0.00500 25.28 (a) These determinations were made in equilibrium with saturated 
Benzene 40 2 0.273 —0.0060 0.00537 28.09 sodium chloride solution whereas the remaining ones were made in 
Benzene 40 3 0.297 —0.0070 0.00584 32.77 equilibrium with dehydrite. 
F-5 60 4 0.058 —0.0012 0.007 19 28.89 
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Fic. 2. Capillary rise of some of the Fulweiler creosotes in 
paper strips of approximately equal benzene calibrations. 


flow for the various papers might be quite similar 
since Whatman No. 4 filter paper was used 
throughout this work. The fallacy of such 
reasoning is immediately evident from the times 
of rise of the benzene calibrations as given in 
Table I. This has made it impracticable to com- 
pare the times of rise of the Fulweiler creosotes 
directly. However, it was possible to select 
several strips whose calibration values were quite 
similar (193+-3 seconds) for a height of rise 
extending from 3 to 10 cm. The corresponding 
height-time data for the penetration of the 
Fulweiler creosotes into these papers have been 
plotted in Fig. 2. These curves which should 
be directly comparable indicate that the rate of 
penetration in paper increases markedly with 
temperatures, the lower boiling creosotes are the 
better penetrants and the differences in the 
rates of penetration among the creosotes are less 
at the higher temperatures. 


THEORETICAL RELATIONS 


Theoretical and empirical considerations of 
experimental data offer a means of approaching 
the problems involved in the capillary penetra- 
tion of fibrous materials with liquids. From a 
scientific viewpoint the applicability of es- 
tablished theoretical relationships to experi- 
mental data should always be considered before 
applying empirical relations. The latter should 
be used only to correlate data which do not 
conform to any of the known exact relations, to 
correlate data when there is no theoretical 
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Fic. 3. dh/dt vs. 1/h for some of the Fulweiler creosotes in 
filter paper. 


equation known or to establish a more simple 
correlation which would be more suited to the 
analysis of a great amount of experimental data. 
This portion of the discussion has therefore been 
limited to a consideration of the applicability of 
established equations of flow in capillaries 
based upon theoretical considerations as ex- 
pressed by Eqs. (3) and (4). 

It seemed more advisable, for reasons which 
become evident later in the discussion, to con- 
sider the applicability of the more complex 
Eq. (4) before Eq. (3). In so doing the method 
suggested by Peek and McLean® has been fol- 
lowed throughout and the values of dh/dt so 
obtained have been plotted against 1/h for some 
of the creosotes as shown in Fig. 3. The resulting 
straight lines are in good agreement with the 
theory developed by the above authors* and 
this agreement was found for all of the creosotes 
studied. 

In Table II are listed the values of slope and 
intercept obtained with each of the creosotes 
tested, together with the values for their cor- 
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responding benzene calibrations. The correspond- 
ing values of Ay cos @ (or 4y times the slope) 
and B (or 8n/dg times the intercept) have been 
computed according to Eq. (4) for each slope and 
intercept and these values are also included in 
Table II. 

The benzene slope values listed in Table II 
indicate that Peek and McLean* must have been 
dealing with paper strips of very similar struc- 
ture since they were able to use a single value of 
benzene slope for the calculation of the relative 
penetration tensions of the pure liquids which 
they investigated. The values of the benzene 
slopes given here vary by approximately 70 
percent from the lowest value found. This 
stresses the importance of calibrating each strip 
of paper with a liquid such as benzene before 
determining relative values of any other liquid 
by this method. 

The values of B (or 8n/dg times the intercept) 
shown in Table II are in reasonably good agree- 
ment for the same paper strip but they vary 
considerably among the different strips. The 
close agreement of the values of B determined 
with benzene and creosote show that this value 
is independent of the liquid used which is in 
accordance with the requirements of Eq. (4). 

It is of interest to note that the values of 1/h 
when dh/dt is equal to zero afford a measure of 
the height to which a liquid will rise in a given 
paper strip providing the structure of the paper 
is assumed to be the same over this height as over 
the portion measured and that the theory can be 
assumed to extend to the ultimate height reached 
by the liquid. It is seen from the curves of Fig. 3 


TABLE III. Values of the ratio of the penetration tensions of 
the Fulweiler creosotes relative to benzene as determined 
in glass and paper. 








SAMPLE 

No. 40°C 60°C 80°C 
1.04( 1.00) 
1.07(1.04) 
0.97(0.94) 


1.05(1.06) 


100°C 
0.97(0.91) 





1.12(1.04) 
1.09( 1.04) 
0.99(0.95) 
1.03(1.03) 
1.04(1.01) 
1.02(1.08) 
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Note: The values of the above ratio which are shown in brackets were 
determined by Eq. (11) which will be discussed later. The remaining 
values were determined by Eq. (4). 

(a) These determinations were made in equilibrium with saturated 
sodium chloride solution whereas the remaining ones were made in 
equilibrium with dehydrite. 
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Fic. 4. Penetration of some of the Fulweiler creosotes into 
paper strips according to Washburn’s equation. 


that the heights to which the creosotes for which 
data are plotted will rise vary from approxi- 
mately 30 to 100 cm. 

The penetration tension values of the Fulweiler 
creosotes relative to benzene have been com-- 
puted from the values of Ay cos @ listed in Table 
II. The ratios of these values to the correspond- 
ing surface tension values relative to benzene are 
given in Table III. In the same table are also 
included values of the same ratio determined by 
an empirical equation but these need not be 
considered at present. The values of the ratio 
listed in Table II indicate that the ratio should 
be equal to one. The variations from unity are 
believed to be due to the limited accuracy of the 
paper strip method. This seems to be indicated 
by the lack of any definite check or trend when 
any one of the creosotes is considered at the 
various temperatures. It is concluded from this 
that the angles of contact for all of the Fulweiler 
creosotes are equal at all of the temperatures 
studied. As it would be very unlikely that all 
of the creosotes would have the same contact 
angle unless the latter were zero, it was also con- 
cluded that their angle of contact against paper 
is either zero or very small. This would also 
follow at once if the former assumption of a 
zero contact angle for benzene against paper is 
correct. 

The plot of #? against ¢ given in Fig. 4 shows 
that the data contained herein do not conform 
to the straight line relationship postulated by 
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Eq. (3). This slight but real curvature may have 
been caused by the small value of h when ¢ was 
zero but it can be demonstrated that the possible 
error incurred in ¢ would not account for the 
curvature found. Washburn obtained this equa- 
tion by considering the hydrostatic head to be 
negligible when dealing with the capillary 
penetration of liquids into very small vertical 
capillaries. It has already been shown that the 
data contained herein conform to Eq. (4) which 
differs from Washburn’s equation only in that 
the hydrostatic head and variations in pore size 
are considered. Since Eq. (4) integrates to give a 
form similar to (3) when hydrostatic head is 
neglected, it appears that the failure of Wash- 
burn’s equation to fit the data must be due to 
neglect of hydrostatic effects. Whitehead and 
Greenfield’ successfully applied Washburn’s 
equation to their data probably because they 
considered only small heights of rise but it must 
follow that their values of r are somewhat in 
error because of neglect of hydrostatic head. 
Moreover, their values of r would not be those 
postulated by Washburn’s equation but would 
be multiplied by the summation of some dis- 
tribution function of 7; i.e., their r would be the 
A of Eq. (4) except for the difference involved 
in not considering hydrostatic head. 

While Eq. (4) is applicable to these data and 
is useful for fundamental research, it is some- 
what difficult to apply in more practical prob- 
lems where many data must be analyzed. It 
seemed advisable because of this to consider the 


applicability of empirical relations to the data in 
hand. 


EMPIRICAL RELATIONS 


The logarithmic diagram for a few of the 
height-time data, as given in Fig. 5, is typical of 
all the experimental results. It is evident from 
these curves that the data obey the empirical 
relation of Bell and Cameron? as given by Eq. 
(1).* The values of m satisfying this equation 
were determined graphically and were found to 


* The latter is based on the assumption that fh is zero 
when ¢ is zero which was not the case for these data 
although the equation was found applicable for the small 
values of h when ¢ was zero. It was found impossible to 
remedy this condition because of the high initial speed of 
penetration especially in the case of benzene. 
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vary from 2.00 to 2.38. These values for m are 
in good agreement with those found by Lucas.‘ 
Moreover, it was found that the values of 
for the creosotes did not differ substantially 
from those of their corresponding benzene 
calibrations. Hence, it must follow that the 
differences in the values of m among the various 
paper strips used must have been caused by 
differences in the structural configuration of the 
paper strips. 

Simmonds® has used the above equation for 
penetration studies but has assumed by analogy 
with Washburn’s equation that k=ry/2n. It 
can be shown from dimensional considerations 
that such an analogy is not justified. 

The similarity of » for the penetration of 
creosote and benzene into the same paper strip 
shows that these liquids obey the same general 
laws in paper capillaries. It should therefore 
follow from Eq. (1) that 


t(benzene) /t(creosote) = C 


(7) 


when the two liquids have been measured over 
the same area of a paper strip and the heights of 
rise are equal. The values of C should be constant 
for each creosote and temperature depending 
upon the physical properties of creosote in paper 
but should be independent of the original size, 
shape and distribution of the capillary pores in 
the paper. The values of C should therefore 
afford a simple and direct measure of the rela- 
tive penetrability of the creosotes in paper 
providing variations among the paper strips are 
only of a structural nature. The latter seems 
reasonable when it is considered that all of the 
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Fic. 5. Logarithmic h vs. t diagram for some of the Fulweiler 
creosotes in paper. 
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strips were not only prepared from the same lot 
of Whatman No. 4 filter paper but also received 
similar pretreatment. 

It is also evident that any pore size changes 
brought about by either liquid would be cor- 
respondingly reproduced in C. The effects of 
benzene on the pore size should be relatively 
constant and need therefore not be considered, 
but the effects of creosotes on pore size may not 
be constant. The latter may be considered by 
some to give erroneous values of C but it is 
believed that this is not the case. In explanation, 
consider two creosotes which differ only in their 
ability to vary pore size. It is evident that the 
creosote maintaining the larger pore size during 
impregnation would penetrate the more rapidly 
and should therefore have the greater penetrat- 
ing power. The values of C should therefore 
portray accurately any differences in penetra- 
bility among the Fulweiler creosotes since original 
structural differences among the paper strips 
have been-eliminated by the benzene calibrations. 
The values of C have been computed and are 
listed in Table IV. 


TABLE IV. Penetrability of the Fulweiler creosotes relative to 
benzene as given by C. 











0.2088 
8 0.1969(a) 








(a) These determinations were made in equilibrium with saturated 
radium chloride solution whereas the remaining ones were made in 
equilibrium with dehydrite. 


The relative penetrabilities of the Fulweiler 
creosotes increase with volatility and tempera- 
ture. A plot of the logio C against the concentra- 
tion of residue distilling above 355°C is shown 
in Fig. 6. The resulting curves are remarkable 
in that they show an inverse proportionality 
between log:.) C at constant temperature and the 
volatility of these creosotes as expressed by 


logie C=(E/c)+F (8) 
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Fic. 6. Relative penetrability—concentration diagram for 
the penetration of the Fulweiler creosotes into paper. 


where c is the concentration of residue not dis- 
tilling up to 355°C and E, F are constants de- 
pending on the temperature. The effects of- 
volatility upon the relative penetrability of the 
Fulweiler creosotes is demonstrated but the 
curves also show the desirability of using high 
impregnation temperatures with these creosotes. 
This is evident since the penetrability not only 
increases with temperature but the differences 
among the creosotes become less. 

In another paper! it was shown that the logy 
viscosity at constant temperature is directly 
proportional to the residue not distilling up to 
355°C in these same creosotes. It should there- 
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. 7. C vs. y/n for the penetration of the Fulweiler 
creosotes into paper strips. 
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TABLE V. Values of G(Cn/y) X 104. 
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SAMPLE 
No. 40°C 


1.806 
1.801 
1.966 
1.820 
1.850 
1.748 
1.969 
1.826 
1.941 
1.830 
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fore follow that the relative penetrability must 
depend upon some inverse function of the vis- 
cosity at constant temperature. Moreover, 
theoretical considerations of liquid flow in 
capillary systems would indicate that the rela- 
tive penetrability should depend upon some 
direct function of the surface tension. In Fig. 7, 
C has been plotted against y/n in accordance with 
the simplest relation suggested by the above 
reasoning. The points shown follow a straight 
line relationship remarkably well and indicate 
that 


C=Gy/n (9) 


where G is a constant for all of the Fulweiler 
creosotes and is independent of temperature. 
The values of this constant are given in Table V. 

The values of G taken as a whole vary quite 
considerably but it is suggested that it would 
be more suitable to consider average values in 
dealing with data which are known to be subject 
to fluctuating errors. When this is done the 
applicability of Eq. (9) to the data is demon- 
strated. 

The relative penetrabilities of the Fulweiler 
creosotes have been shown to depend upon their 
respective viscosities and surface tensions in 
accordance with Eq. (9). It must follow accord- 
ingly that the times of ‘rise of the Fulweiler 
creosotes into the same or identical paper strips 
should be directly proportional to their vis- 
cosities and inversely proportional to their 
surface tensions. Hence, since it has been shown 
that h™ is proportional to ¢ when 7/y is held 
constant and ¢ is proportional to »/y when h is 
held constant, it must follow that 


h"=Kyt/n, (10) 
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where 7 should be a constant for each liquid and 
paper depending upon the hydrostatic head as 
will be shown later and K should be a constant 
for each paper strip depending upon the struc- 
tural configuration of the paper. It is evident 
from theoretical considerations that the de- 
pendence of K upon the structural configuration 
of the paper strips is based on the assumption of 
a zero or constant paper-liquid contact angle. 
Although this is apparently the case for the 
Fulweiler creosotes as was shown previously, 
it would evidently not always be true. Since the 
contact angle is known to reduce merely the 
effects of surface tension during penetration, the 
more general equation of flow may be written 

h=K’'y cos 6t/n, (11) 
where K’ would be dependent only upon the 
structural configuration of the paper. 

An equation relating the factors affecting the 
penetration of the Fulweiler creosotes into paper 
strips has been derived from purely empirical 
considerations of experimental data. The ratios 
of the surface tensions of these creosotes de- 
termined in glass to the corresponding values of 
y cos 6 determined in paper by the application of 
Eq. (11) are given in Table III. The values are in 
good agreement with those determined by Eq. (4) 
and the same conclusions of zero or small contact 
angle creosote-paper can be drawn from the 
results. It is also evident that the use of Eq. (11) 
for the determination of penetration constants 
is much the simpler of the two methods. 

It will be noted that if Eq. (11) is considered 
in a differential form it will differ from Eq. (4) 
essentially in that no hydrostatic head term 
appears and that the / in the denominator of 
Eq. (4) has become /”~' in the differential form 
of Eq. (11). It is concluded from this that the 
value of m empirically corrects Eq. (11) for the 
hydrostatic head term. 


MOISTURE EFFECTS 


The effects of different amounts of adsorbed 
water on the penetration constants of the Ful- 
weiler creosotes were only considered at two 
moisture contents. It was felt that the close 
agreement in the values obtained did not war- 
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rant a complete investigation of the effects of 
water below fiber saturation. Moreover, un- 
published data on other liquids have shown that 
varying amounts of adsorbed water do not 
change the values of the penetration constants. 
It must be remembered in this connection that 
any changes in pore size caused by the adsorbed 
water are eliminated from consideration since 
the height-time data for both the benzene 
calibration and creosote were determined on 
papers having the same moisture content. 

It is believed that differences in the penetra- 
tion constants of the Fulweiler creosotes might 
be found in the presence of free water for in this 
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case interfacial tension effects would be intro- 
duced. The interfacial tensions of these creosotes 
have been shown! to vary greatly so that the 
contact angles between creosote and _ paper 
containing free water should also vary. The 
penetration tension (y cos #) would vary accord- 
ingly and consequently the rate of penetration 
into paper strips. The inclusion of a study of the 
rate of penetration into paper strips containing 
free water was not possible since methods for the 
duplication of conditions are not available. This 
is essential since the pore size of the capillaries 
must be taken into consideration by calibration 
with another liquid. 
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Heat Generation and Anisotropy of Rubber Compounds in the Flexometer 
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Data are given showing that it is possible for rubber to 
soften in one direction and stiffen in another direction 
simultaneously during a test in the St. Joe Flexometer. Both 
effects are mainly the result of heat, the softening being 
due to plastic flow, while the stiffening is due to the sub- 
stantial elimination of plastic flow because of the rapidity 
of flexure. Static and dynamic tests on heated test pieces, 
and stress-strain data on various sections of test specimens 
at several stages of breakdown indicate that failure in 
dynamic tests is due partially to mechanical fatigue, but 


mainly to heat effects. Results on anisotropic pigments such 
as asbestine, clay, magnesium carbonate, and acicular and 
non-acicular zinc oxides show that the heat generation, 
breakdown time, flexing load, and vertical deflection are 
functions of the shape and orientation of the pigment 
particles. Variations in breakdown time from 16 to 160 
minutes have been obtained by changing the orientation 
of pigment particles. Data on ‘‘grain”’ effect in high gum 
stocks are also given. 





INTRODUCTION 


N a previous paper! a new machine for evaluat- 

ing the breakdown characteristics of rubber 
compounds when flexed under compressive 
loads was described. This machine, which is 
known as the “St. Joe Flexometer,” differs from 
previous flexometers in that the flexing forces 
can be measured at all times, and the end point 
is definitely indicated on the machine. With this 
machine it was shown that when rubber is 
flexed at a constant horizontal deflection under 
either constant vertical compression or constant 
vertical load an increase in flexing force is re- 
quired to maintain the deflection at a constant 
value. No attempt was made in the previous 
paper to explain the phenomenon of the simul- 
taneous occurrence of plastic flow and stiffening 
of rubber test specimens which was indicated by 
the data obtained. It was simply pointed out 
that the data showed good correlation with ac- 
tual road tests on commercial solid tire and 
carcass stocks. While a general relationship 
between heat generation and breakdown time 
has been indicated! ? a survey of the literature 
did not disclose any information regarding the 
phenomena taking place in rubber when flexed 
under compressive loads. In the present paper a 


‘ Havenhill and MacBride, Ind. Eng. Chem., Anal. Ed. 
7, 60 (1935). 


* Barnett and Mathews, Ind. Eng. Chem. 26, 1292 
(1934). 
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study has been made of the simultaneous oc- 
currence of plastic flow and stiffening action 
taking place, which has logically led to the 
investigation of anisotropy. 


DESCRIPTION OF ST. JOE FLEXOMETER 


Fig. 1 shows a view of the entire machine. 
Fig. 2 is a close-up view of the machine showing 
the sample under vertical compression and 
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Fic. 1. View of the entire machine. 
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The horizontal load, vertical deflection or com- 
pression, and temperature at the center of the 
test specimen have been plotted against time, 
For these and subsequent tests, unless otherwise 
noted, the following formula was used: 


Parts by Weight 













Rubber 100.0 
Sulphur 4.0 
D.O.T.G. 1.4 
Gas Black 9.0 
Zinc Oxide * 150.0 


Cure 65 min. at 141.1°C 


* Electrothermic process zinc oxide of medium particle size was used. 
Average diameter’ by number =0.26 micron 
Specific surface = 1.83 m?2/g 





























It will be noted that as the test progresses the 
Fic. 2. Close-up showing a sample under vertical com- Vertical compression increases (height of test 
pression and horizontal deflection. specimen decreases), which means that the 

sample is softening and flattening out. The 


horizontal deflection. Fig. 3 is a drawing of the horizontal load or flexing load, however, in- 
machine in which the various parts have been ‘'©#5©S- This means that _ the test ite eae 
laiacled the force to produce a given deflection in the 
The test specimen, which is a cylindrical plug 
3.81 cm in diameter and 3.81 cm in height, is 
: placed between the two parallel Bakelite face z £4 
plates and under compression by means of the = ©.5 
vertical load weight 6, acting through lever arm : 3 2 af 
7. The upper face plate is then rotated and the 5 — S 
sample, which acts as a coupling, causes the C — 5 
lower plate to rotate. The lower face plate is — <—=—(28) 
mounted on the horizontal roller bearing car- S : 
riage traveling on the track and can be moved off 
center from the upper plate any given amount | : —(4) 
by manipulation of handwheel 20. The load @ J &\, — 6 
required to produce or maintain this horizontal I3, 
deflection is read on scale 19. The horizontal A Scale Dial 
deflection itself is read on dial 26. By a slight 19 
rotation of handwheel 20 the linkage 16 can be 
lengthened or shortened, thus decreasing or 
increasing the horizontal load and _ horizontal 


© 2, 20 














Fic. 3. Diagram of machine. 








° “3 ae 1. Rubber test specimen 15. Horizontal load weights 
deflection on the sample. In this manner it is 3. Upper face plate Sh, ‘Slates Gee elatite 
possible to run tests under constant horizontal 3. Lower face plate 17. Bell crank 

P ° 4. Centering jig 18. Link 
deflection (HD) or under constant flexing load 5. Hydraulic jack ae ‘aici hail 
(HL). For a more complete description of the 6. a load weights 20. Handwheel for deflection ad- 

° o« ° 7. Lever justment 
machine refer to the original article. 8. Upper vertical shaft 21. Parallel linkage 
Fig. 4 shows what happens during a constant 9. V-belt pulley , 22. Double divided scale 
h e 1 d fl ° Th - 1 1 d . 10. Lower vertical shaft 23. Straight edge for micrometer 
orizonta efiection test. e vertical loa 1S 11. Carriage 24. Pointer for vertical compression 
maintained constant at 250 kg and the horizontal 12. Carriage rollers 25. Pointer for horizontal deflection 

° ° . ° 13. Track 26. Dial-type micrometer gauge 

deflection is maintained constant at 0.673 cm. 44. Hydraulic jack 37. Canons ee 
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Fic. 4. Constant horizontal deflection test. Vertical load, 
250 kg; horizontal deflection, 0.673 cm. 
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horizontal plane increases, indicating a marked 
stiffening of the stock. The temperature also 
increases rapidly during the test, which is 
largely due to the combined effect of mechanical 
working and resultant frictional heat; the low 
heat conductivity and diffusivity of the rubber 
compound; and, in part, to the Joule effect. 
Previous tests' carried out under constant 
vertical compression and constant horizontal 
deflection conditions show that the increase in 
horizontal load is not due to a change in the 
shape of the specimen, but is due primarily to a 
structural change or stiffening of the rubber. 
Evidence of this is presented in Fig. 4A. 


. 


SIMULTANEOUS STIFFENING AND PLASTIC FLOW 
OF RUBBER 


To show how it is possible for the rubber to 
both soften and stiffen in different directions at 
the same time during a flexometer test, it was 
felt that if stress-strain data could be obtained 
for various sections of plugs run for different 
lengths of time in the flexometer, the data would 


OF RUBBER COMPOUNDS 


60 





s 





3 








8 
Vertical Load in Kg 





Horizontal Load in Kg 


s 


























Ki 
Time in Minutes 


Fic. 4A. Constant vertical compression test. 
Vertical compression, 1.61 cm. 
Horizontal deflection, 0.714 cm. 


throw some light on what was happening during 
test. The difficulties encountered in slicing, 
buffing, and dieing out tensile strips from the 
regular test specimens or plugs led to the use of a 
laminated plug. It was found that a specimen 
built up of disks (3.81 cm diameter and 0.19 cm 
in thickness) died out from vulcanized sheets 
would act as a unit in the flexometer, giving 
results which checked the results obtained by 
using the regular molded cylindrical plug. 

Figs. 5, 6 and 7, show tensile strength, 100 per- 
cent modulus (stress at 100 percent elongation), 
and elongation at break data for laminated 
plugs run to 1/10, 1/2 and 9/10 breakdown in the 
flexometer. The tensile strength of the stock 
before flexing was 185 kg/cm*; modulus at 100 
percent elongation, 46 kg/cm*; and elongation 
at break, 380 percent. The results show that 
during the first part of the test there is a slight 
increase in tensile strength toward the center of 
the plug. This is brought about by the cold 
working of the stock at relatively low tempera- 
ture and is probably due to the formation of a 
new grain structure. As the test progresses the 
tensile strength gradually falls off due to the 
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Fic. 5. Tensile strength at 1/10 breakdown time. 
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Fic. 6. Tensile strength at 1/2 breakdown time. 


combined effect of heat, cure and mechanical 
working. These effects are greatest toward the 
center of the plug. 

While these data show a slight increase in 
tensile strength during the first part of the flex- 
ometer test, they do not explain the increased 
stiffening of the specimens as the test progresses. 
Since the stiffness (as shown by the increase in 
horizontal load) is measured in the flexometer at 
high temperature, and the tensile strength was 
measured at low temperature, it was felt that the 
increased stiffening might be due to the effect 
of heat. In order to check this idea, plugs were 
heated for 20 minutes and then tested in the 
flexometer. 

In Fig. 8 the horizontal load has been plotted 
against the time of test for samples preheated at 
107°C, and tested while hot and after cooling. 
The curve for a specimen not preheated is in- 
cluded to show that the additional cure due to 
preheating causes a lower horizontal load. Similar 
curves for specimens preheated at other tem- 
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Fic. 7. Tensile strength at 9/10 breakdown time. 
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Fic. 8. Constant horizontal deflection tests on preheated 
specimens. Vertical load, 250 kg; horizontal deflection, 
0.610 cm. 


peratures also indicate that the higher flexing 
load is a heat effect rather than a cure effect, 
and is a function of the temperature of preheat- 
ing. In Fig. 9 the temperature of preheating has 
been plotted against the time of breakdown 
only, and it will be noted that as the tempera- 
ture increases the breakdown time decreases. 
This relationship holds for tests in which the 
rate of heat generation is much greater than the 
rate of heat dissipation, as is usual test procedure 
in the flexometer. In the case of stocks (under- 
cured stocks especially) tested under very low 
vertical loads and very low horizontal deflections 
so that the rate of heat generation approaches 
that of dissipation, this relationship may not 
hold true. These results indicate that the stiffen- 
ing of the rubber is mainly a heat effect and that 
heat has a predominant influence on failure 
during flexing. 

Since the heat generated by mechanical 
working and possibly the Joule effect has caused 
a stiffening in the horizontal plane, a similar 
stiffening effect might be expected in the vertical 
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Fic. 9. Temperature preheating vs. breakdown time. 
Vertical load, 250 kg; horizontal deflection, 0.610 cm. 
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Fic. 10. Static compression test on hot and cold plugs. 
Vertical load, 250 kg; horizontal deflection, 0.714 cm. 


plane. However, the results of flexometer tests 
show that during the test the samples actually 
soften and flatten out, and assume a permanent 
set. Static compression tests on heated plugs 
show that for the first ten seconds or so there is 
a slightly increased stiffness due to heat, but 
that plastic flow soon sets in and causes the 
flattening out and resultant permanent set. 

Fig. 10 shows the data of such a static com- 
pression test run on a hot and cold plug. The 


OF RUBBER COMPOUNDS 














4, Permanent Set 


























st $—d 


Time in Minutes 


Fic. 11. Permanent set acquired during flexometer test. 
Vertical load, 250 kg; horizontal deflection, 0.673 cm. 


first plug was heated 20 minutes at 157°C, 
cooled to room temperature, and then tested in 
the flexometer without rotation under 250-kg 
vertical load and 0.714-cm horizontal deflection. 
This plug showed slight evidence of plastic 
flow but there was no measurable permanent 
set at the end of 20 minutes. The other plug was 
heated 20 minutes and immediately placed in the 
flexometer under the same static conditions. 
The curves show that at the beginning of the 
test of the hot sample there is a stiffening in both 
the vertical and horizontal planes, as indicated 
by the initial high horizontal load and low 
vertical compression as compared with the cold 
sample. After the first few seconds of test 
plastic flow sets in as shown by the rapid in- 
crease in vertical compression (flattening out of 
the sample) and decrease in horizontal load. This 
plastic flow results in a permanent set. Fig. 11 
shows permanent set data plotted against time 
of test. These data were obtained by stopping 
several dynamic tests at different stages of 
breakdown. Since it was found that the set 
measured after a 24-hour rest period was very 
nearly the same as that measured on the hot 
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plugs immediately after the test, only the set 
after 24 hours is shown. In the dynamic or 
regular flexometer test the specimen is flexed so 
rapidly in the horizontal plane (875 r.p.m.) 


that plastic flow in that plane cannot take place 


to any substantial degree. In the vertical plane 
the compressive force acts continuously in the 
same direction, so that flow does take place. 
This rapid flexing and reversal of stress in the 
horizontal plane without flexing in the vertical 
plane makes it possible for the particular rubber 
compounds tested to exhibit the simultaneous 
stiffening and plastic flow during the flexometer 
test. 

These data show very definitely that rubber 
is stiffer at high temperatures than at low tem- 
perature when plastic flow is substantially 
eliminated. From thermodynamic considerations 
also*: * we would expect the rubber to be stiffer 
at increased temperatures. While ordinary stress- 
strain data® usually show a lower modulus at 
high temperatures, stress-strain tests made at 
very high speeds® show a higher modulus for 
temperatures up to 140°C. Higher modulus 
could undoubtedly have been obtained on the 


layers of laminated plugs (Figs. 5, 6 and 7) if the 
stress-strain tests had been made at a higher 
speed than is customary with the standard stress- 
strain machine, and at the temperatures reached 
in the flexometer test. 


ANISOTROPY 


Since rubber compounds when tested in the 
flexometer exhibit simultaneous plastic flow and 
stiffening in different directions, it might be 
expected that pigments having different particle 
shapes would, when incorporated into the rubber 
stocks, produce anisotropy of the compounds and 
exert a pronounced effect on heat generation. 
Since no mention is made in the literature regard- 
ing the effect of anisotropy on heat generation, 
the following data are presented. 

It was first thought that tests could be made 
on cubes built up from layers each having the 


3 Gerke, Ind. Eng. Chem. 22, 73 (1930). 

we and Snyder, Rubber Chem. Tech. 8, 151 
(1935). 

5 oo and Cope, India Rubber World 79, 64 
(1928). 

6 Williams, Ind. Eng. Chem. 21, 872 (1929). 
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Fic. 12. Test on cubic plug. Vertical load, 290 kg; horizontal 
deflection, 0.760 cm. 


grain in one direction. These cubes could then 
be tested with any one of the three axes in the 
vertical plane of the machine to show the differ- 
ences in breakdown characteristics due to 
“grain” effects and orientation of the pigments. 

It might be mentioned here that the resistance 
to shear of a block of wood along the grain is less 
than that across the grain, and much less than 
that perpendicular to the grain. The horizontal 
load might be considered analogous to the shear- 
ing force, and the facing of the test piece in the 
flexometer would determine whether shear should 
occur with, across, or perpendicular to the 
“grain.” 

Fig. 12 shows results of tests on a 3.81 cm cube 
of the same stock as that used for tests shown in 
Fig. 4. Numbers at the end of the curves indicate 
the top face of the specimen during test. Because 
of unequal stresses in the different directions, 
considerable vibration was encountered and 
duplicate results were very hard to obtain. The 
cubic test specimen was then discarded and the 
molded cylindrical specimen was used. To align 
the pigment particles in different directions, three 
methods of preparing the samples prior to 
vulcanization were used. Method 1: The uncured 
stock is sheeted out after remilling for two 
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Fic. 13. Pictorial representation of the orientation of a 
rectangular plate-shaped pigment particle. 


minutes on a laboratory mill (roll temperature 
49°C) to 0.381 cm thickness. After a 2-hour rest 
period, strips 5.08 cm wide, with a length equal 
to the width of the sheeted stock are cut off 
across the grain and rolled up, making the grain 
of the rubber parallel to the axis of the cylinder. 
Method 2: In this method strips are cut in the 
direction of the grain of the sheeted stock and 
rolled into cylinders. Here the pigment particles 
are in a spiral perpendicular to the axis of the 
specimen. Method 3: In this method disks died 
out of the uncured stock are built into a lami- 
nated specimen having the pigment particles in 
layers perpendicular to the axis. 

Fig. 13 is an exaggerated pictorial representa- 
tion of the orientation of a rectangular plate- 
shaped pigment particle subjected to the three 
methods of preparing samples. 

Microscopic examination showed that pig- 
ment particles were lined up during milling with 
their longest dimension in the direction of 
milling and their shortest dimension normal to 
the rolls. Examination of the stocks to be tested 
showed that this alignment was present after 
vulcanization. If we consider a pigment having 
needle- or rod-shaped particles, these particles 
would be aligned, in method 1, with the long 
dimension parallel to the axis of the cylinder. 
In the flexometer the needles would be vertical, 
and would offer greatest resistance to the hori- 
zontal load or flexing force. During the constant 
horizontal deflection test, we would expect a 
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high horizontal load and a relatively short 
breakdown time due to the rapid heat generation. 
In methods 2 and 3 the small dimension of the 
particles would be parallel to the axis of the test 
specimen, offering a low resistance to the 
horizontal load or flexing force, and we would 
consequently expect lower flexing loads and 
longer running time in the flexometer test. 

If we consider plate-shaped particles with 
square or round faces, we would expect high 
flexing loads with short breakdown times in 
methods 1 and 2 where the larger dimension of 
the particle is parallel to the vertical axis, 
offering greater resistance to the flexing load. 
In method 3 the plates would form a laminated 
structure in which the thickness of the particle 
would be parallel to the vertical axis, and we 
would expect a much lower horizontal load and 
much longer breakdown time. 

If the pigment particles are rectangular 
plates, we would expect method 2 to give a lower 
flexing load and longer breakdown time than 
method 1. Method 3 would give much lower 
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Fic. 14. Constant deflection tests on stock containing 
clay. Vertical load, 250 kg; horizontal deflection, 0.570 cm. 
Curves 1, 2 and 3 indicate method of preparation. 
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Fic. 15. Constant deflection tests on stock containing 
asbestine. Vertical load, 250 kg; horizontal deflection, 0.661 
cm. Curves 1, 2 and 3 indicate method of preparation. 
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Fic. 16. Constant deflection tests on stock containing 
acicular ZnO. Vertical load, 250 kg; horizontal deflection, 
0.635 cm. Curves 1, 2 and 3 indicate method of preparation. 
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Fic. 17. Constant deflection tests on stock containing 
nonacicular ZnO. Vertical load, 250 kg; horizontal deflec- 
tion, 0.661 cm. Curves 1, 2 and 3 indicate method of 
preparation. 


flexing load and considerably longer running 
time than method 2. 

In order to check this theory, stocks were 
made up by using the formula referred to above 
with all but one volume of zinc oxide replaced 
by equivalent volume loadings of magnesium 
carbonate, clay, asbestine, acicular zinc oxide and 
nonacicular zinc oxide. These stocks were tested 
in the flexometer with a constant vertical load 
of 250 kg and constant horizontal deflections 
indicated on the graphs. 

Figs. 14, 15, 16, 17 and 18 show the horizontal 
load-time curves obtained in the tests. The 


temperature at the center of the test specimen 
after 10 minutes of flexing is indicated along 
the curves. The method of preparation of the 
plug, which indicates the grain direction, is given 
by the number at the end of each curve. A 
high-gum stock (the above formula without 
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Fic. 18. Constant deflection tests of stock containing 
magnesium carbonate. Vertical load, 250 kg; horizontal 
deflection, 0.508 cm. Curves 1, 2 and 3 indicate method of 
preparation. 
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Fic. 19. Constant deflection tests on pure gum stock. * 
Vertical load, 215 kg; horizontal deflection, 1.130 cm. 
Curves 1, 2 and 3 indicate method of preparation. 


pigment) is shown in Fig. 19 for the purpose of 
comparison. Temperature measurements were 
not taken during tests of this stock in order to 
eliminate any possible tendency of the hypo- 
dermic-needle thermocouple to cause premature 
failure. Samples of high-gum stocks generally 
have a tendency toward a mechanical type of 
failure. Table I includes time, vertical compres- 
sion, horizontal load, and temperature data for 
each test at the start, after 10 minutes, and at 
the breakdown point. 

Results of the flexometer tests on stocks 
containing clay (Fig. 14) show that the break- 
down characteristics of test specimens prepared 
by methods 1 and 2 are very similar, while in 
those prepared by method 3 the horizontal load 
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is lower and the running time considerably 
longer, with slow heat generation. These data 
indicate that the particles are thin plates. 

In the case of stocks containing asbestine 
(Fig. 15) the differences in breakdown charac- 
teristics between specimens prepared by methods 
1 and 2 are very pronounced. This indicates that 
the length of the particle is considerably greater 
than the width. The very low flexing load, long 
time, and slow heat generation of specimens 
made by method 3 indicate that the thickness of 
the particle is considerably less than either the 
length or width. In effect, the particles of as- 
bestine act as if they were long thin needle-like 
crystals or ribbons. 

The curves for stocks containing acicular 
zinc oxide (Fig. 16; the oxide used being an 
American process product having an average 


diameter’ of 0.36 micron and a specific surface 
of 1.54 m?/g) show that specimens made by 
method 2 run for a much longer time under a 
lower flexing load than specimens made by 
method 1. These data indicate that the length 
of the particle is much greater than the width, 
and the similarity of the curves for specimens 
made by methods 2 and 3 show that the thick- 
ness is about the same as the width, indicating 
a close approach to the true needle shape. 

For nonacicular zinc oxide (Fig. 17) the data 
indicate a more nearly isotropic shape. 

The results of these tests enable the prediction 
of particle shapes which check fairly closely 
those shown by microscopic examinations of the 
pigments used. However, breakdown character- 
istics of plugs made by method 3 show greater 


7 Perrott and Kinney, J. Am. Ceramic Soc. 6, 417 (1923). 


TaBLeE I. Test of stock with various pigments. VC =vertical compression. HL = horizontal load. Values are given at starting 
time, at the end of 10 minutes and at the time of breakdown. 
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; METHOD 
ERI F PE : 
or a TIME VC HL Temp. — vuapa TIME VC HL Temp. 
PIGMENT TEST RATION (min.) (cm) (kg) (°C) PIGMENT TEST RATION (min.) (cm) (kg) (°C) 
1 0 1.16 7.0 30 1 0 1.16 8.0 30 
Start 2 0 1.14 5.5 30 Start 2 0 1.02 7.5 30 
3 0 1.19 0.0 30 3 0 1.09 5.0 30 
ia 1 10.0 1.38 14.8 142 NONACICULAR 1 10.0 1.13 12.7 118 
HD =0.570 cm 10-min. 2 10.0 1.32 11.5 132 ZINC OXIDE 10-min. 2 10.0 1.04 12.2 110 
—— 3 10.0 1.23 1.0 85 HD =0.661 cm 3 10.0 1.14 7.8 102 
1 27.5 2.14 38.0 216 1 320 1.70 360 177 
Break 2 345 2.12 3328 213 Break 2 41.0 1.62 33.0 176 
3 141.0 2.03 23.0 205 3 66.5 1.70 315 175 
1 0 1.08 7.0 30 1 0 1.40 9.0 30 
Start 2 0 1.04 6.0 30 Start 2 0 1.35 7.0 30 
3 0 1.08 3.0 30 3 0 1.32 0.0 30 
ASBESTINE : 1 10.0 1 35 15.0 168 MAGNESIUM 1 10.0 1.94 25.1 1 74 
HD =0.661 cm 10-min. 2 10.0 1.05 8.0 110 CARBONATE 10-min. 2 10.0 1.86 21.8 168 
? 3 10.0 1.10 3.0 88 HD =0.508 cm 3 10.0 1.52 1.0 121 
1 17.0 1.75 24.0 213 1 16.0 2.24 30.5 208 
Break 2 41.0 1.82 20.0 218 Break 2 16.5 2.21 29.5 207 
S 157.5 1.85 220 196 3 40.0 2.20 21.2 192 
1 0 1.14 5.5 30 1 0 1.75 —8.0 30 
Start 2 0 1.04 3.4 30 Start 2 0 1.71 —6.5 30 
3 0 1.12 2.0 30 3 0 1.68 —7.0 30 
ACICULAR 1 10.0 1.50 21.8 129 —— 1 10.0 1.68 —5.7 — 
Zinc OXIDE 10-min. 2 10.0 1.14 5.5 104 HD =1.130 10-min. 2 10.0 165 —-40 — 
HD =0.635 cm 3 100 1.17 34 94 at iain aes 3 100 163 -45 — 
16.5 1.93 46.3 200 1 96.0 1.98 265 — 
Break 2 45.3 1.90 32.0 198 Break 2 1065 1.909 235 — 
61.5 1.98 38.2 182 3 159.0 1.90 230 — 
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differences from characteristics of methods 1 and 
2 than would be predicted by the theory and 
microscopic studies. These differences would 
indicate. anisotropy of the rubber hydrocarbon 
itself, and the curves for the high-gum stock 
(Fig. 19) show evidence of such anisotropy which 
may possibly vary in its effect with different 
pigments in compounded stocks. This long run- 
ning time of specimens prepared by method 3 is 
also in part due to the fact that the stock is 
approaching the point at which heat dissipation 
is as rapid as heat generation. 

As a further check on the theory, flexometer 
tests were run on a stock containing magnesium 
carbonate. From the curves of Fig. 18, a thin 
plate-shaped particle was predicted. This shape 
differs from the usual one which is needle-like, 
but a microscopic examination revealed that the 
plate type of particle was actually predominant. 

While the data in general check the theory, 
which is based on particle shape, there is evi- 
dence of other factors which enter into the picture 
so that quantitative predictions cannot be made. 
Several of these factors are ‘“‘skin”’ effect in the 
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laminated plugs, possible lack of complete 
alignment of particles in the milling operation, 
and variable anisotropy of the rubber itself in 
different compounds. There is also microscopic 
evidence of a breaking up of some of the pigment 
particles during test. 


CONCLUSIONS 


Data have been presented to show that it is 
possible for rubber compounds when tested in the 
flexometer to exhibit simultaneously a plastic 
flow in one direction, and a marked stiffening 
in another direction resulting from the elimina- 
tion of plastic flow by rapid flexing. Both effects 
are shown to be due mainly to heat generation 
or temperature effects. The rate of heat genera- 
tion is in part dependent on the anisotropy of the 
rubber compound, which is much greater than is 
indicated by some of the ordinary physical tests 
such as tensile, tear and abrasion. The theory 
that the degree of anisotropy is influenced by the 
shape and orientation of the pigment particles 
and of the anisotropy of the rubber itself has 
been in a measure experimentally verified. 
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The stresses around large obstacles in strained rubber were studied with the assumption that 
the effects would be similar with actual pigment particles. The effects of variations in the size 
and shape of the obstacles on the stress concentration were investigated and an explanation 
offered for the stiffiening of rubber by pigments based on the change in slope of the rubber 
stress strain curve. The various average diameters of a series of zinc oxide samples were meas- 
ured and their relationship with some physical properties of rubber compounds determined. 





te importance of flow in the rubber on the 

reinforcing properties of pigmented systems 
has been emphasized by Park.! He suggests: (1) 
that in the presence of a finely divided pigment, 
the flow which occurs when a piece of rubber is 
stretched, takes place in the capillary spaces be- 
tween the pigment particles; (2) that some modi- 
fication of the laws of liquid flow may govern the 
behavior of rubber with reference to pigments 
embedded in it; and (3) that the forces causing 
increased stiffness in pigmented rubber com- 
pounds are similar to those causing increased 
resistance to flow of liquids in tubes of capillary 
dimensions. Thus increasing fineness of subdivi- 
sion and the resulting fineness of capillary spaces 
between the particles should be accompanied by 
an increase in reinforcing properties. 

It would be desirable to study the actual 
stresses around pigment particles in rubber under 
strain but so far no suitable microscopic set-up 
has been devised. A few years ago the writer, 
resorting to analogies, measured the strains and 
stresses around large particles with the assump- 
tion that the strains would be relatively the same 
with small particles. For this study holes of the 
desired size and shape were cut in strips of 
calendered but uncured gum stock and fitted 
with pieces of an uncured semihard rubber com- 
pound. After vulcanization squares were marked 
on the tensile sheets as shown in Fig. 1. 

With no particles in the sheet the strains re- 
mained uniform on elongation as shown by Fig. 2. 
However, when particles were inserted, as in 
Figs. 3 to 7, the stresses became uneven and 





* Investigator, Research Division, The New Jersey Zinc 
Company, Palmerton, Pa. 
1C. R. Park, News Ed., Ind. Eng. Chem. 11, 345 (1933). 


strains were greater around the particles. These 
uneven strains caused by the particles resemble 
very much the lines of flow of a liquid past an 
object. The percentage elongation was obtained 
from the photographs and, with the load known, 
it was possible to draw the stress strain curve for 
any part of the tensile sheet. 

In Fig. 8, stress is shown for varying distances 
from the particle for circular and rectangular 
particles. When no particle is present the cor- 
responding stress in the rubber is 125 pounds per 
square inch, corresponding to 225 percent elonga- 
tion, and this stress is shown by the horizontal 
line. All of the curves approach this line as the 
distance from the particle increases. 

Two conclusions are possible from Fig. 8. 


1. Stresses are greater at the surface of large particles 
than smaller ones. This checks Schippel? who found that for 
the same material the volume increase under strain varied 
directly with the particle size. 

2. The stress in the vicinity of a rectangular particle 
is great if the particle is oriented with its length parallel 
to the direction of pull but small if oriented normal to this 
direction. 


The stress strain curve of unpigmented rubber 
may be considered as of two parts, one parallel 
to the strain axis (soft) and the other parallel 
to the stress axis (stiff). According to Figs. 3 to 8, 
the addition of pigment causes the strain on con- 
siderable of the rubber to be raised to the stiff 
part of the stress strain curve when the average 
rubber would normally be in the soft part of the 
curve. 

The magnitude of this stiffening action may 
be determined by the size, shape, size distribution 
or dispersion of the pigment. Fig. 9 shows a series 


2 Schippel, Ind. Eng. Chem. 12, 33 (1920). 
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Fic. 1. Marking of vulcanized squares. 





of stress strain curves for circular and rectangular 
particles. The values for the stress are those 
found by extrapolating curves similar to those in 
Fig. 8 to the stress axis. The elongations, as 
abscissae, are those for the rubber as a whole. 
Here again the greatest stiffening is obtained 
with the rectangular particle oriented lengthwise 
to the direction of extension. In the case of the 
circular particles, the magnitude of the stiffening 
by the particle depends on its diameter. Between 
particles of one and two centimeters in diameter, 
there is an increase in strain of 32 percent for the 
section of maximum concentration of stress; 
between particles of one and three centimeters in 
diameter, the increase is 48 percent. These in- 
creases are in the ratio of one to one and one-half 
which is the same as the ratio of diameters. It 
has been calculated* that for an expanded tetra- 
hedral arrangement (which gives a greater inter- 
particulate distance than an expanded cubical 
arrangement) in a 20-volume carbon black stock 
the particles are distant from each other slightly 
more than their own radius. Yet Figs. 3 to 7 
have shown that the field of influence of the 
particles extends well beyond the distance repre- 
sented by the radius of the particles; thus these 
effects would tend to be relatively greater in a 
rubber compound containing a moderate amount 
of pigment. 

From the information accumulated in work 
with actual pigments in rubber, two general 
theories have developed. The one based on 
abrasion, tear and stress-strain studies assumes 
a heterogeneous structure to be responsible for 


* Wiegand, India Rubber J. 73, 31 (1927). 


Fic. 2. Uniform strain, no particles in the sheet. 


reinforcement by increasing stiffness. An attempt 
to give the mechanism of this action has been 
made in this paper. The second based largely on 
stress strain measurements and heat of swelling 
data explains reinforcement as a surface phenom- 
enon. Among those attacking the problem 
from the first angle Haslam‘ has endeavored to 
determine which of the various average diameters 
of a pigment was most significant in determining 
the physical properties of rubber and how great 
an error would result if an average diameter 
other than the correct one were used. Mixtures of 
four zinc oxide samples of different particle size 
were prepared so that several samples had the 
same. surface area but differed in number of 
particles per unit weight while others differed in 
surface area but had the same number of particles 
per unit weight. Stress strain properties, abrasion 
and tear resistance were measured on the rubber 


* Haslam, The Rubber Age (N. Y.), Nov. 10, 1932. 
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Fics. 3 to 7. Nonuniform strains produced by the insertion of particles. 


compounds. In general, as the specific surface or 
the number of particles were increased the rein- 
forcing properties improved with the better 
relationship in the case of specific surface. How- 
ever, even in the samples in which specific surface 
was held constant, the physical properties varied 
somewhat and these variations were ascribed to 
variations in uniformity since they were propor- 
tional to the variations between the particle 
sizes of the component oxides. Table I gives the 
various diameters, the specific surface, number 
of particles per gram and uniformity coefficients 
of a series of pigments and pigment mixtures 
prepared according to the method of Haslam as 
well as the equations used in their calculation. 
A lengthy discussion of these diameters is not 
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Fic. 8. Stress in pounds per square inch vs. distance from 
edge of particle. 


justified here since the subject has been covered 
by Perrot and Kinney’ and by Green.* d; is the 
arithmetic mean diameter, d; the diameter from 
which specific surface may be calculated, D the 
diameter from which the number of particles per 
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Fic. 9. Stress vs. strain for various shapes of particles. 


5 Perrot and Kinney, J. Am. Ceramic Soc. 6, 417 (1923). 
® Green, J. Frank. Inst. 204, 713 (1927). 
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TABLE I. Various diameters, the specific surface, number of particles per gram and uniformity coefficients of a series of pigments 
and pigment mixtures. Equations defining these quantities in terms of the density (), the number of particles 
of various sizes (m), and the fraction weights of components of mixtures (w). 

















= = _ Ind _ Xdyw _ (=nd\+ ,_ (zw/d 4 =nd* 
Ss pd;' N pD3’ dy =n’ Ud,= zw ’ a=( 2) : v= (=) . 3a 
. =w Ind\i .. =w \} =nd‘ Udw 
Uds =w/d;’ - (= ’ uD=(= 7) ' d= snd Ua= lw 
SAMPLE MIXTURE OF Ss N X1010 dy Ud A U d3 Uds D UD ds Ud 
A ) 5.08 3975 0.132 1.39 0.143 1.55 0.21 1.48 0.165 1.64 0.27 1.25 
B single 2.36 448 .256 1.16 314 1.42 45 1.29 .342 1.38 53 1.20 
ec pigments 1.57 107 44 1.36 .50 1.45 .675 1.4 55 1.50 17 1.345 
D 55 4.5 1.28 1.40 1.44 1.53 1.92 1.49 1.59 1.65 2.22 1.385 
E A+B 3.64 2013 .147 1.12 .173 1.21 .298 1.3 .207 1.36 .36 1.54 
F A+C 3.64 2420 .138 1.15 .160 we | .293 1.56 .196 1.43 476 1.31 
G A+D 3.64 2720 .133 1.28 .150 1.43 .293 1.51 .187 1.60 .89 .985 
H A+C 3.30 2014 .140 1.09 .165 1.22 aan 1.35 .207 1.37 52 1.215 
I A+D 2.84 2013 .133 1.18 .154 1.33 .374 1.32 .207 1.51 1.23 .90 
M B+D 1.57 255 .187 1.01 .226 1.10 .392 1.21 .272 1.23 50 1.2 
N A+B+D 1.57 606 .162 .87 .208 1.00 .467 1.14 .272 1.14 .67 1.11 
L? A+D 1.57 890 .136 91 .172 1.06 .683 .98 .292 1.23 1.79 74 
J A+B 2.70 890 .270 1.02 .323 .955 .693 1.35 .292 .93 1.27 1.11 
K A+C 2.27 890 161 .84 .209 .985 .680 1.08 .309 1.13 1.55 .875 








gram may be calculated, A is the diameter of the 
particle of average surface and d, the diameter of 
the particle of average weight. The uniformity 
coefficients U, were calculated from the following 
equation given by Green: 


U=d\(¥in/2¥-n(d,—d)?)}, 


in which d; is any average diameter and is intro- 
duced to make U independent of magnification 
and dependent upon the relative values of parti- 
cle size and the sum of the squares of the “‘resid- 
uals,’’ 2(d,;—d)*. When all the particles are equal 
in size, U becomes infinite and the average 
diameters coincide. For mixtures the equation is 
expressed as a function of the fractional weight 


(W): 


0-a(e5/2|SfEs2a-]) 


d is used to designate any fixed diameter of the 
mixture. 

An inspection of Table I shows that the d,; and 
A diameters have little significance in this 
instance, as for example sample J which is a 
50.6 : 49.4 mixture of sample A (d;= : 132; 
A=0.143) and sample D (d,=1.28; A=1.44) has 
a d; diameter of 0.133 and a A diameter of 0.154. 
These diameters are actually smaller for sample 
I than for sample H which is a 49.3 : 50.7 mixture 


of samples A and C. The reason is that these 
diameters are so markedly affected by the fine 
particles and not appreciably by the coarse. 
Further, the values of the uniformity coefficients 
of the various mixtures are not in agreement 
with the uniformity which would be expected 
from the particle sizes of the components of the 
mixture. The agreement is satisfactory when the 
coefficient is calculated on the d, diameter and 
the next best results are obtained when the d; 
diameter is used, yet even in this case samples 
G(A+D), and F(A+C) appear more uniform 
than sample E(A+B). 

Lack of space makes it impossible to present 
all the data obtained on the properties of these 
pigments in rubber compounds but Figs. 10, 11 
and 12 show some of the results obtained in a 
rather unusual test which measures principally 
the flow properties of the compounds. The form- 
ula used was as follows: 


Smoked sheet 100 
Sulphur 4 
Diphenylguanidine 2 
Zinc oxide 30 Vol. 


These results were obtained with a pendulum 
which has been described elsewhere’? and which 
measures: (a) resilience, or the percentage of the 
impact energy returned to the pendulum by the 


' 7Barnett and Mathews, Ind. Eng. Chem. 26, 1292 
(1934). 
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Fic. 10. Resilience and indentation vs. number of particles 
per gram. 
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Fic. 11. Resilience and indentation vs. specific surface 
of particles. 


rubber and (b) indentation, the depth of pene- 
tration of the hammer-shaped head of the pen- 
dulum into the sample. Figs. 10, 11 and 12 show 
the results for resilience and indentation plotted 
against the number of particles per gram, the 
specific surface and the d, diameter, respectively ; 
the optimum cure for resilience was used in each 
case. 

The agreement between the particle size 
measurements and the physical properties is 
obviously better in the last two cases than in the 
first and the agreement with the d, diameter 
seems somewhat better than with specific sur- 
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Fic. 12. Resilience and indentation vs. diameter (d,) 
of particles. 
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Fic. 13. Abrasion index vs. diameters d3; and d,. 


face. The spread between samples having the 
same d, diameter is as great as between those 
having the same surface but in the specific sur- 
face curve there are definitely three groups of 
samples which if viewed independently would 
result in the conclusion that resilience and inden- 
tation increased with increasing specific surface, 
a conclusion not warranted in view of the data as 
a whole. , 

In Fig. 13, the data given by Haslam com- 
paring the abrasion index with specific surface 
have been recalculated to show a comparison of 
this property with the d; and d, diameters. The 
d, diameter seems to give a somewhat better 
agreement with the physical property than does 
the d; diameter or the specific surface calculated 
from it. For instance, the series of pigments at 
constant surfaces show a definite trend in size as 
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indicated by the d, diameter which is in line with 
the abrasion index. 

Perrot and Kinney in their discussion of the 
various average diameters have given a partial 
explanation of these results. They believe that 
the d; average diameter is the logical one to 
apply to paint pigments and fillers since the 
uniformity of the materials is so high that very 
little error is introduced. ‘‘Consider, however, the 
pigment ‘ranging from ultramicroscopic particles 
to grains as coarse as sand.’ The specific surface 
of such a pigment would be considerably less 
than the specific surface indicated by the average 
diameter, because the relative surfaces of the 
particles per unit weight of material have been 
neglected in calculation of this value.’’ It seems 
very probable that the range of sizes in these 
synthetic mixtures of pigments involves this 
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error. A number of investigators have endeavored 
to determine the effect of blending pigments. 
Dawson* reports that neither combinations of 
clay and carbon black nor of two zinc oxides 
follow the mixture law in tensile properties, or 
hysteresis, either in extension or rebound tests. 
Lefcaditis and Cotton® found a departure from 
the mixture law but in the opposite direction 
while Beaver and McKay" who also worked with 
clay and carbon black obtained values in con- 
cordance with the mixture law. If it were possible 
to measure the proper average diameters of these 
pigments in the state of dispersion existing in the 
rubber compounds, it might be possible to 


explain some of these apparent discrepancies. 


8 Dawson, Trans. I.R.I. 2, 92 (1926). 
® Lefcaditis and Cotton, Trans. I.R.I. 8, 364 (1932). 
© Beaver and McKay, Ind. Eng. Chem. 23, 294 (1931), 
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Force at a Point in the Interior of a Semi-Infinite Solid 


RayMonpD D. MINDLIN, Department of Civil Engineering, Columbia University 
(Received March 13, 1936) 


A solution of the three-dimensional elasticity equations for a homogeneous isotropic solid 
is given for the case of a concentrated force acting in the interior of a semi-infinite solid. This 
represents the fundamental solution having a singular point in a solid bounded by a plane. 
From it may be derived, by a known method of synthesis, the solutions for the semi-infinite 
solid which correspond to the solutions known as nuclei of strain in the solid of indefinite 


extent. 





1. INTRODUCTION 


NE of the fundamental results in the theory 

of elasticity is the Kelvin solution! for a 
force applied at a point in a solid of indefinite 
extent. It is well known that stress distributions 
for a number of problems of practical importance 
may be obtained from the Kelvin solution by 
methods of synthesis and superposition.” For 
example, the problem, first solved by Lamé, of 
the stresses in a spherical container under uni- 
form internal or external pressure may be 
treated by superposing stresses derived from 
Kelvin’s results. First we form, from the Kelvin 
single force solution, the stresses for a double 
force, that is, a pair of equal and opposite forces 
acting at neighboring points. By superposing 
three mutually perpendicular double forces, we 
arrive at the solution known as the center of 
compression. Finally, the addition of a uniform 
tension or pressure to the center of compression 
leads to the solution of Lamé’s problem. 

The classical problem of Boussinesq dealing 
with a normal force applied at the plane boundary 
of a semi-infinite solid* has found practical ap- 
plication in the study of the distribution of 
foundation pressures, contact stresses, and in 
certain problems of soil mechanics. It is known 
that Boussinesq’s problem nfay also be solved 
by superposing solutions derived from Kelvin’s 
results. To do this we form a line of centers of 
compression by integrating the stresses for a 
single center of compression. By combining this 


1Sir W. Thomson, Cambridge and Dublin Math. J. 
(1848). See also A. E. H. Love, Mathematical Theory of 
Elasticity, fourth edition (Cambridge, 1927), p. 183; and 
ne, Theory of Elasticity (New York, 1934), p. 

? Love, p. 186, and Timoshenko, p. 323. 

* J. Boussinesq, Applications des Potentials . . . (Paris, 
1885). Also Love, p. 191, and Timoshenko, p. 331. 


result with the original Kelvin single-force solu- 
tion, the Boussinesq formulas are obtained. 

The Kelvin solution may be used in studying 
stresses due to a force applied at a great distance 
from a boundary while the Boussinesgq solution is 
applicable in the case where the force acts at the 
surface. The solutions described in this paper fill 
in the gap between the two by giving the stresses 
for the case where the force is applied near the 
surface. Such a condition is approached in a 
number of practical problems such, for example, 
as in the case of a guy wire anchor. Again, by 
integrating these new solutions along a line, we 
may approximate the conditions produced by a 
friction pile or an anchor rod. 

Solutions of the elasticity equations such as the 
single force, double force, center of compression, 
and line of centers of compression are referred to 
as nuclei of strain.2 The number of nuclei of 
strain which may be obtained by synthesis from 
the Kelvin solution is unlimited and further 
combinations lead to the solutions for such 
important cases as Southwell’s problem of the 
spherical ¢avity in an unlimited solid under 
simple tension* and Cerruti’s problem of a force 
applied tangentially at the plane boundary of a 
semi-infinite solid.§ 

Each of the nuclei of strain obtained from 
Kelvin’s results is a solution applicable to a solid 
of indefinite extent. With the solution for a 
force in the interior of a semi-infinite solid, as 
described in this paper, we may apply the method 
of synthesis to develop a new series of nuclei of 
strain for which the stresses vanish on the plane 
boundary of the semi-infinite solid. Such solutions 
may be termed half-space nuclei of strain. 


926) V. Southwell and H. J. Gough, Phil. Mag. 1, 71 
1926). 

5 V. Cerruti, Acc. Lincei. Mem. fis. mat., Roma 13, 81 
(1882). Also Love, p. 241. 
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The notion of such half-space nuclei of strain 
suggests some interesting possibilities. By em- 
ploying these half-space nuclei as unit solutions, 
it may be possible to find solutions analogous to 
those obtained from the Kelvin nuclei but with 
the additional feature that a new plane boundary 
is introduced. For example, the solution anal- 
ogous to Southwell’s problem would pertain to a 
spherical cavity near the surface of a semi- 
infinite solid; the solution analogous to the 
problem of Boussinesq would be that for a force 
applied at the surface of a body bounded by two 
perpendicular planes (i.e., the transmission of 
force through ‘‘quarter-space’’); the solution 
analogous to the hollow sphere under uniform 
pressure would apply to an eccentrically hollow 
sphere. 

The present paper deals with the fundamental 
solution for the single force in the interior of a 
semi-infinite solid. The problem is divided into 
two parts: (1) force normal to the boundary; 
| (2) force parallel to the boundary. 

The discovery of the solutions followed a 
study of H. M. Westergaard’s interpretation® 
of the problems of Kelvin, Boussinesq and 
Cerruti in terms of Galerkin vectors and the 
results obtained are shown to include these 
problems for certain limiting conditions. 





2. THE GALERKIN VECTOR 


The components of stress in an elastic iso- 
tropic solid were expressed by B. Galerkin’ in 
terms of partial differential coefficients of three 
functions satisfying a fourth-order equation. 
P. F. Papkovitch® and H. M. Westergaard® have 
expressed the displacements and stresses in 
terms of a vector stress function whose scalar 
coefficients are the Galerkin functions. The dis- 
placements are given, in Westergaard’s form of 
the Galerkin vector, by® 


u=(1/2G)[2(1—pw)AX —(0/dx) div F], 
v=(1/2G)[2(1—pn)AY—(0d/dy) div F], 
w= (1/2G)[2(1—")AZ—(0/dz) div F], 


*H. M. Westergaard, Lectures on Elasticity while 
Visiting Professor at the University of Michigan (Summer 
Session, 1934). 

7B. Galerkin, Comptes rendus 190, 1047 (1930). 

(1932) F. Papkovitch, Comptes rendus 195, 513, 754 
1932). 
' 935) M. Westergaard, Bull. Am. Math. Soc. 41, 695 
1935). 


(1) 
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where: X, Y, Z are the Galerkin functions, 
i, j, k constitute an orthogonal system of unit 
vectors 
F=iX+jY+kZ is the Galerkin vector, 
A is Laplace’s operator, 
G is the modulus of rigidity, 
p is Poisson’s ratio. 


The stresses are given by® 
o,=2(1—p)dAX/dx+ (uA — d2/ dx?) div F, 
o,=2(1—p)dAY/dy+ (ud — d2/ dy?) div F, 
o,=2(1—y)dAZ/dz+ (uA — */dz") div F, 

Tyz=(1—y) (dA Y/dz+ 0AZ/dy) 

— (d*/dydz) div F, 

T2z=(1—p)(0AZ/dx+ 0AX / dz) 

— (0?/dzdx) div F, 

try =(1—p)(GAX /dy+ dA Y/dx) 

— (d?/dxdy) div F. 


(2) 





For the case of zero body force, the Galerkin 

vector must satisfy the biharmonic equation, 
A’F =0. (3) 

It was observed by Professor Westergaard‘ 
that the coefficient of a single axially symmetrical 
component of the Galerkin vector is identical 
with Love’s function'® for strain symmetrical 
about the corresponding axis. 


3. NUCLEI OF STRAIN IN TERMS OF GALERKIN 
VECTORS 


The solutions given in sections 4 and 5 were 
obtained by superposition of nuclei of strain for 
an unlimited solid. It has been found con- 
venient, in applying the method of superposition, 
to express the nuclei of strain in terms of Galerkin 
vectors. 

In the following list of functions, A, B, 
od - H are constants and R=(x?+y?+2*)!. 
It may be observed that each of the functions 
satisfies the biharmonic equation and that the 
stresses derived from them vanish at infinity. 


10 Love, p. 274; Timoshenko, p. 309. 











FORCE IN THE 
GALERKIN VECTORS FOR NUCLEI OF STRAIN IN AN 
UNLIMITED SOLID 
A. Single force (Kelvin’s problem). 
(1) iAR (single force in x direction) 
(2) jAR (single force in y direction) 
| (3) KAR (single force in z direction) 
B. Double force 
(1) iBx/R (double force in x direction) 
(2) jBy/R (double force in y direction) 
(3) kBz/R (double force in z direction) 
C. Double force with moment 
(1) iCz/R (double force in x direction with 
about y axis) 
(2) jCx/R (double force in y direction with 
about 2 axis) 
(3) kCy/R (double force in z direction 
about x axis) 
(4) iCy/R (double force in x direction 
about 2z axis) 
(5) jCz/R (double force in y direction with moment 
about x axis) 
(6) kCx/R (double force in z direction 
about y axis) 


moment 
moment 
with moment 


with moment 


with moment 


D. Line of double forces with moment 

(1) iDz log (R+-x) (type C(1) along x axis from x 
tox=— «) 

(2) jDx log (R+y) (type C(2) along y axis from y 
to y=— «) 

(3) kDy log (R+<2) (type C(3) along z axis from z 
to z=— «) 

(4) iDy log (R+x) (type C(4) along x axis from x 
to x=— ») 

(5) jDz log (R+y) (type C(5) along y axis from y 
to y=— «) 

(6) kDx log (R+2) (type C(6) along z axis from z 
to z=— 0) 


ll 
r—) 


ll 
o 


ll 
o 


ll 
o 


ll 
o 


ll 
Oo 


E. Center of compression or dilatation 
(1) i£ log (R+x) 
(2) jE log (R+~y) (All three yield the same stresses) 
| (3) KE log (R+2) 
F. Line of compression or dilatation of constant strength 
(1) iF[x log (R+x)—R] (along x axis from x=0 to 
x==— 2 ) 
(2) jFLy log (R+y)—R] (along y axis from y=0 to 
y= - 2 ) 
(3) kF[z log (R+2)—R] (along z axis from 2z=0 to 
s=— 0) 
G. Doublet (double center of compression-dilatation) 
(1) iG(1/R) (axis of doublet parallel to x axis) 
(2) jG(1/R) (axis of doublet parallel to y axis) 
(3) kG(1/R) (axis of doublet parallel to z axis) 
H. Linearly varying line of doublets with strength propor- 
tional to distance from the origin 
(1) iH[z log (R+z2)—R] (type G(1) along z axis from 


z=0 to z=— ~) 
(2) jH[x log (R+x)—R] (type G(2) along x axis from 
x=0 to x=— ) 
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(3) kH[y log (R+y)—R] (type G(3) along y axis from 
y=0 to y= — ~&) 

*(4) iH[Ly log (R+y)—R] (type G(1) along y axis from 
y=0 to y= — ~) 

(S) jH[s log (R+s)—R] (type G(2) along s axis from 


z=0to s=— ~) 
(6) kH[x log (R+x)—R] (type G(3) along x axis from 
x=0tox=— 0) 


4. ForcE NORMAL TO THE BOUNDARY OF A 


SEMI-INFINITE SOLID 


The semi-infinite solid is considered to be 
bounded by the plane z=0, the positive z axis 
penetrating into the body. A force P is applied 
at point (0, 0, +c) and acts in the positive z 
direction (Fig. 1). 

Since the stress will be symmetrical about the 
2 axis, we need consider only the k component of 
F. Transforming to cylindrical coordinates 
r, 0, z and noting that Z is independent of 6, we 
find, from Eqs. (1) and (2), the following ex- 
pressions for the displacements and stresses in 
symmetrical cylindrical coordinates :'° 


U = —(1/2G)(a°Z/araz), 

w= (1/2G)[2(1—w)AZ — 02Z/ 027], 
o,= (0/02) AZ — 0°Z/ dr], 

o9= (0/02) AZ —(1/r)(dZ/dr)], 
o:=(0/02)[(2—p)AZ— 072/027], 
trz=(0/dr)[(1 —w)AZ — 0°Z/ 027]. 


Phone z+0 

















Fic. 1. Force normal. to the boundary in the interior of a 
semi-infinite solid. 
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We first apply the Kelvin solution at point 
(0, 0, +c) by replacing R with R,=(r?+(z—c)?)! 
in solution A(3) of section 3. We note that 
stresses o, and r,, are produced on the plane 
z=0. The boundary conditions for z=0 are 


Loz J:mo=[7rz]2-0=0. (5) 


It is therefore necessary to add further solutions 
to remove the boundary stresses and these solu- 
tions must introduce no new singularities in the 
region z=0. Nuclei of strain of types A(3), B(3), 
E(3), F(3) and G(3), applied at (0, 0, —c), are 


RAYMOND D. 








} 


MINDLIN 


found to satisfy these conditions if the proper 
values are assigned to the constants A, B, E, F 
and G. 

To transfer nuclei of strain from the origin to 
(0, 0, —c), we replace z with (z+c) and R with 
R2=(r+(z+c)?)! in the formulas of section 3, 

The values of the constants are obtained from 
the boundary Eqs. (5) and the equilibrium 
condition 


--f 2nro.dr, (z>C). (6) 
0 


This procedure leads us to the required Galerkin vector: 
F =[Pk/8r(1—y) ]{Rit(3—4y)R2—2c(s+c)/R2—4(1—2y)c log (Re+2+c) 


+4(1—)(1—2n)[(s+c) log (Ret+s+c)—ReJ+2c?/Re}. (7, 


The solution for a force in the interior of a semi-infinite solid and normal to the boundary is seen 
to be compounded of six nuclei of strain which, in an unlimited solid, represent: (1) a single force at 
(0, 0, +c); (2) a single force at (0, 0, —c); (3) a double force at (0, 0, —c); (4) a center of compression 
at (0, 0, —c); (5) a line of centers of compression extending from z= —c to z= — ~ ; (6) a doublet at 





(0, 0, —c). 
Collecting terms in Eq. (7) we obtain :"! 


F =[Pk/87(1—y)]{|Rit+[8u(1 —v) —1]R.+4(1—2u)[(1—w)s—pe] log (Re+z+c)—2cz/Re\. (8) 


Eqs. (8) and (4) yield the displacements and stresses: 


















































; Pr rs—c (3—4y)(2-—c) 4(1—p)(1—2u) 6c2(s+c) 
"TerG(—w Re’ ROR Rits+0 ORS | 
P r3—4u 8(1—p)?—(3—4yu) (2-—c)? (3—4y)(e+c)?-—2c2 6ce2(2+c)? 
~ 16"G(1—wk Ri Ry R:3 Rs TRS } 
P (1—2yu)(2—c) (1—2u)(s+7c) 4(1—p)(1—2yu) 3r?(ze—-—c) 
————— Rs! R(Rets+0) RY 
6c(1 —2yu)(s+c)?*—6c7(s +c) —3(3—4y)r2(z—c) 30cr?s(s+c) P 
ocniciace eet 
P(i—2p)f(s—c) (3—4y)(s+c) —6c 4(1—p) 6c(s+c)* 6c7(s+c) 
anal Ri R3 R(Ret+s+c) Rs <a 
P (1—2y)(s—c) (1—2p)(s—c) 3(z—c)? 
¢->=- _| - - a ~ 
8x(1—yp) Ri} R;' Ri 
3(3—4y)2(s+c)?*—3c(s+c)(52—c) 30c2(2+c) 
_ Rs R# } 
Pr 1—2u 1-2 3(s-—c)? 3(3—4y)2(s+c)—3c(32+c) 30c2(s+c)? 
eer Re RY ORS Rs R¥ } 


" R. D. Mindlin, Comptes rendus 201, 536 (1935). 
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Fic. 2. Normal stress on planes parallel to boundary forthe Fic. 3. Vertical deflection of the surface for the case of a 
case of a force normal to the boundary (u=0.3). force normal to the boundary (u=0.3). 


Fig. 2 shows the distribution of normal stress on horizontal planes, and Fig. 3 shows the vertical 
deflection of the surface of the semi-infinite solid. 

When c—~ all terms in Eqs. (9) containing R, vanish and the solution becomes that for Kelvin’s 
problem where the force is applied at (0, 0, +c) in the positive z direction. The corresponding Galerkin 
vector 1s 


F=k[ PR,/8r(1—y) ]. (10) 


When c—0 Eqs. (9) give the stresses and displacements for Boussinesq’s problem.' These stresses 
are represented by the Galerkin vector 


F=k(P/2r)[2uR+(1—2u)s log (R+2)]]. (11) 


The solution for the problem of Boussinesq was given by Professor Westergaard® in the form of the 
Galerkin vector 


F=k(PR/2r) (12a) 
combined with the potential function 
= —[(1—2u)/2m]P log (R+3). (12b) 
The relation between ¢ and F is given by 
grad @¢=2(1—yw)AF—grad div F. (13) 


As a further verification of the solution, we may compare the stresses with those obtained for the 
corresponding two-dimensional case by E. Melan.™ While the approach to the solution in three 
dimensions was simplified through the use of cylindrical coordinates, the comparison with Melan’s 
formulas is facilitated by referring the components of stress to rectangular coordinates. Eqs. (9) then 


become: 





"The stresses and displacements for the problems of 13 E. Melan, Zeits. f. angew. Math. und Mech. 12, 343 
Boussinesq and Cerruti are given by F. Vogt, Det Norske (1932). This comparison was suggested by Professor S. 
egy ae t Oslo, Math.-Naturv. klasse, No. 2, Timoshenko. 

Pp. 9, ; 
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P. [~ —2yu)(z—c) 3x*(z—c) Pu. —2u)[3(s—c) —4y(z+c) } 

































































'eei—wl ORR R 
3(3—4y)x?(z—c) —6c(2+c)[(1—2p)z—2yc] 30cx?s(z+c) 
Rg RY 
a 
R2(Re+2+¢) Rx(Ro+2+c) R-? 
P “(1—2y)(z—c) 3y%(z-—c) (1—2y)[3(s—c) —4y(s+c) ] 
ee R® Re | Rs 
3(3 —4y) y?(s—c) —6c(2+c)[(1—2p)s—2ye] 30cy*s(s+c) 
OO R:5 7 Ri 
4(1—y)(1—2y) y? y? 
Ro(Re+2+0) “gaa ™ 
- SS ee 
8x(1—p) R;3 R;' ORY 
3(3—4y)2(2-+c)?—3c(e+c)(5z—c) 30cz(z+c) 
a Rs RS | 
Py (1—2u) (1-—2u) 3(s—c)? 3(3—4y)2(2+c)—3c(32+c) 30c2(2+c)*7 
aor Re RP ORS Rs Re 
Px - (1-—2u) (1-—2y) 3(s-—c)? 3(3—4y)2(2+c)—3c(32+c) 30c2(z+c)*] 
S(i-oL Re Re ORY R: Re OS 
Pxy [ 3(s—c) 3(3—4y)(s—c) 4(1—n)(1—2y) 1 1 30c2(z+c) 
ww Mik Re RS R2?(Ro+2+c) rare _ R,' } 


We consider a uniform distribution of forces 
of magnitude p per unit of length along a line 
through (0, 0, +c) parallel to the y axis. If db is 
a small element of this line at a distance } from 
the z axis, the stresses due to the uniform pressure 
p acting on db are found by substituting pdd for 
P and (y—d) for y in Eqs. (14). If, now, we 
integrate these formulas with respect to )} 
between the limits — © and + ~, we obtain the 
corresponding solution for plane strain. By taking 
account of the relations between the elastic 
constants for plane strain and plane stress and 
also noting an interchange of axes, this pro- 
cedure leads directly to the stresses given by 
Melan. 





5. FoRCE PARALLEL TO THE BOUNDARY OF A 
SEMI-INFINITE SOLID 


In this case there is no axial symmetry and we 
employ rectangular coordinates (x, y, z), and 
Eqs. (1) and (2). The force P is applied at (0, 0, 
+c) and acts in the positive x direction (Fig. 4). 

We again begin by applying the Kelvin solu- 
tion at (0, 0, +c), this time using type A(1). To 
remove the stresses on plane z=0, introduced 
by the single force, we employ nuclei of strain of 
types A(1), C(6), D(6), G(1) and H(1), again 
substituting (2+c) for zand R2= (22+ y?+ (z+c)’)! 
for R. The constants associated with these 
solutions are determined from the boundary 


conditions [o,]-0=[Ty:le-0=[Tez]Jeo=0 (15) 
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and the equilibrium condition 


p=-4f [ Td ydx, 
0 0 


(z>Cc). 


The result of these operations is the Galerkin 


vector! 


F=([P/8r(1 —p) |(if R,+R2—2c?/Re 


+4(1—)(1—2u)[(2+c) log (Re+2z+c)—Re]} 


+k} 2cx/R.+2(1—2u)x log (R2+2+c)}). 


Thus, for the case of a force in the interior of a 
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semi-infinite solid, parallel to the boundary, the 
solution is obtained by superposition of six nuclei 
(16) of strain for the unlimited solid: (1) a single force 
at (0, 0, +c); (2) a single force at (0, 0, —c); 
(3) a doublet at (0, 0, —c); (4) a semi-infinite line 
of doublets extending from z=—c to z=—« 
with strength proportional to the distance from 
z=-—c; (5) a double force with moment at 
(0, 0, —c); (6) a semi-infinite line of double 
forces with moment extending from z=—c to 
ym eo 
The displacements and stresses, obtained from 
Eqs. (1), (2) and (17), are found to be: 
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Fic, 4. Force parallel to the boundary in the interior of a 
semi-infinite solid. 


Fig. 5 shows the distribution of normal stress on 
horizontal planes, and Fig. 6 shows the vertical 
deflection of the surface of the solid along the 
x axis. 

When c—~ all terms in Eq. (18) containing 
Rz vanish and the solution becomes that for 
Kelvin’s problem with the force applied at 
(0, 0, +c) in the positive x direction. The cor- 
responding Galerkin vector is 


F=i|PR,/8r(1—y)}. (19) 


When c—0 Eqs. (18) give the stresses and 
displacements for the problem of Cerruti.” 
These results are represented by the Galerkin 
vector 


F=(P/4r(1—p) }(i{ R+2(1—p)(1 — 2p) 
X([z log (R+s) —R]} 
+k} (1—2yu)x log (R+2)}). (20) 


The solution for the problem of Cerruti was 
given by Professor Westergaard® in the form of 
the Galerkin vector 


F=[P/4x(1—,) ][iR 
+k(1—2y)x log (R+z)] (21a) 
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Fic. 5. Normal stress along lines parallel to the x axis on 
planes parallel to the boundary for the case of a force 
parallel to the boundary (u=0.3). 





Fic. 6. Vertical deflection of the surface along the x axis for 
the case of a force parallel to the boundary (4 =0.3). 


combined with the potential function 
@=[(1—2n)P/2e][x/(R+z)]. (2b) 


Following the same procedure as in the case of 
the normal force, the stresses given in Eqs. (18) 
may be reduced to Melan’s formulas for the two- 
dimensional analog.” 

The author is indebted to Professor H. M. 
Westergaard for valuable criticisms of the 
manuscript. 
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